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ABSTRACT 

1  In  Multi-Scale  Profiler  (MSI*),  a  fm'ly  falling  dropsondc.  has  been  us<*d  over 
tin-  la>>l  12  years  to  resolve  oceanic  slu'ar  variance  at  vertical  scales  from  a  few  hundred 
meters  down  to  nearly  a  centimeter.  Because  MSI*  yielded  the  first  complete  oceanic 
'hear  spectra,  it  is  important  to  document  the  methods  by  which  they  were  produced. 
Large  »<  ales  are  measures!  by  an  electromagnetic  current  meter  ( KCM  ).  microscales  bv 
airfoil  probes,  and  int<TiiM*diate  scab's  by  an  acoustic  current  meter  (  ACM  ).  The  ACM 
detect'  velocity  relative  to  the  instrument,  so  the  platform  motion  must  be  known 
to  determine  the  water  velocity.  Primarily,  tin  ACM  measurements  are  affected  In 
tilt  osc  dial lotis  and  by  the  gross  (point -mass)  motion  of  the  vehicle;  the  former  is 
intern'd  from  accelerometer  data,  and  tin*  latter  is  constructed  from  a  model  of  the 
vehicles  response  to  oc cam  shear.  Horizontal  forcing  on  the  array  of  diag  finishes 
and  turning  blade'  at  the  tail  complicates,  the  response  bv  causing  MSI*  to  react 
strongly  to  line  t ual  ions  with  sc  ah's  near  t In-  Hist rument  's  length  of  1 .1  in  We  examine 
the  effect  sol  thi'  response  on  spec  Ira  of  the1  ACM  measurement',  noting  partuularlv 
a  deep  notch  near  ().‘2  c  pm  (cycles  per  meter).  In  account  for  such  spectral  features, 
the  model  of  Hayes  et  al.  for  the  1  Ol’S  dropsondc  was  modified  so  that  it  correctlv 
patalneteri/ed  eur  large  tail  force.  We  discuss  the  dynamics,  data  processing,  and 
model  formulation  relevant  to  produc  tion  of  oceanic  velocity  profile's  from  the  ACM 
data,  and  present  analytic  transfer  func  tions  derived  from  Lourter  transform'  of  the 
model  equations  which  guide  selection  of  optimal  values  for  the  model  parameters. 
Velocity  profiles  and  shear  spectra  resulting  from  the'  motion-corrected  ACM  data 
compare  well  with  KCM  results  at  large  sc  ales.  Owing  to  MS!*  s  high  sampling 
frequencies,  the  ACM  data  provide  a  ditvet  cheek  on  the'  verae  it  y  of  the  airfoil  probe 
data  at  scale's  of  0. 1  1  rn.  The  overall  strengl  h  of  our  result  s.  as  well  as  t  he-  weaknesses 
and  uncertainties,  is  summarizc'd. 
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I.  INTRODUCTION 

Observations  during  the  past  decade  demonstrated  the  ability  of  the  Multi-Scale 
Profiler  (MSP)  to  resolve  velocity  and  temperature  profiles  in  the  upper  1  km  of 
the  ocean.  The  need  for  MSP  became  apparent  after  analysis  of  the  Fine  and  Mi¬ 
crostructure  Experiment  (FAME)  of  1975,  in  which  different  profilers  resolved  various 
portions  of  the  vertical  wavenumber  spectra  (Gregg  and  Sanford,  1980,  Gargett  et 
al.,  1981).  Although  much  was  learned  during  FAME,  measuring  viscous  dissipation 
rates,  e,  with  one  vehicle,  fine-scale  shear  with  another,  and  diffusive  dissipation  rates, 
\,  with  yet  another  proved  insufficient  for  relating  shear  to  turbulence.  In  response, 
we  began  designing  MSP  in  1980;  we  first  used  it  in  the  open  ocean  in  1983  and  have 
taken  it  on  five  subsequent  cruises  (Table  1).  Observations  from  CSALT,  PATCHEX, 
and  Tropic  Heat  2  yielded  the  first  fully  resolved  vertical  spectra  taken  by  a  single 
instrument  (Gregg  et  al.,  1991,  1993),  and  profiles  in  the  Florida  Strait  took  the 
combined  measurements  close  to  bottom  topography  (Winkel  et  al.,  1992).  Here  we 
report  MSP  characteristics  and  data  processing  affecting  our  published  results.  In 
addition,  because  the  need  for  fully  resolved  spectra  will  continue  past  the  life  of 
MSP,  we  document  key  features  of  the  design. 

To  measure  velocity  at  scales  from  1  kin  to  nearly  10  mm,  MSP  combines  three 
previously  developed  techniques:  electromagnetic  (Sanford  et  al..  1978),  acoustic 
travel  time  (Evans  et  al..  1978:  Hayes  et  al.,  1984),  and  hydrodynamic  lift  on  small 
airfoils  (Osborn.  1974).  By  sensing  voltages  induced  as  sea  water  flows  through  the 
earth's  magnetic  field,  the  electromagnetic  current  meter  (ECM)  yields  velocity  pro¬ 
files  largely  unaffected  by  the  response  of  MSP  to  ocean  currents.  However,  these 
profiles  are  not  absolute  in  that  they  include  unknown,  depth  independent  offsets. 
Proper  treatment  of  the  ECM  measurements  requires  that  MSP  rotate,  and  ECM 
results  are  reliable  only  at  scales  greater  than  the  distance  for  one  such  rotation. 
The  Neil  Brown  acoustic  current  meter  (ACM),  owing  to  its  stability  and  low  noise, 
can  resolve  scales  close  to  the  0.2  tu  path  length  between  its  transducers.  Because 
the  ACM  detects  motion  relative  to  the  vehicle,  minimal  platform  motion  is  desir¬ 
able.  I  nfortunatelv.  the  design  of  MSP  makes  it  highly  responsive  to  fluctuations 
with  scales  close  to  or  greater  than  its  4.3  m  length.  In  converting  the  relative  mea¬ 
surements  to  ocean  velocities,  correcting  for  tlx-  vehicle  response  requires  a  model. 
The  airfoil  or  shear  probes  sense  velocity  fluctuations  relative  to  MSP.  Their  data 
are  transformed  into  shear  spivtra  to  estimate  viscous  dissipation  rates.  Also,  they 
extend  MSP's  resolution  down  to  10  mm:  smaller  scales  are  smoothed  by  the  probe 
tips. 

To  avoid  vehicle  dynamics  associated  with  a  rapid  spin,  we  use  rotation  lengt  hs  of 
15  40  in.  much  larger  than  the  scales  directly  generating  most  turbulence  in  the  ther 
mocline.  Given  the  KCM's  limitation,  we  must  rely  on  the  ACM  and  the  associated 
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Table  1.  MSP  cruises  with  nominal  fall  and  rotation  rates,  Qrot  and  VF0>  from  the 
main  thermocline  (below  2-3  MPa).  In  1990,  the  smaller-diameter  tail  array  replaced 
the  larger  array,  leading  to  faster  fall  rates. 


Cruise 

Name 

Date 

Location 

Rotation  Rate, 
Drot  (Hz) 

Fall  Rate, 

Wo  (ms"1) 

PAC83 

May  1983 

35°  N,  131°W 

Tested  several  configurations 

CSALT 

Nov  1985 

12°N,  56. 5° W 

0.011,  0.018 

0.27 

PATCH EX 

Oc.t  1986 

34°N,  127° W 

42°N,  126°W 

0.017 

0.25 

Tropic  Heat  2 

Apr  1987 

0°N,  140°W 

1,  2,  6,  &  11.5°N 

0.007 

0.002,  0.007 

0.27 

Florida  Straits 

Jun  1990 

27°N,  79.1  >.9°W 

0.013 

0.32 

COARE  3 

Feb  1993 

1.7°S,  0,3°N,  156°E 

0.009 

0.30 

vehicle  response  model  to  resolve  the  scales  bridging  internal  waves  and  turbulence. 
Results  from  the  model  compare  favorably  with  those  from  the  ECM  at  scales  larger 
than  the  rotation  length.  The  smallest  scales  detected  acoustically,  0.1-1  m,  require 
little  correction  and  provide  a  check  on  the  airfoil  signals,  which  can  be  contaminated 
in  this  band  by  low-frequency  noise  and  temperature  sensitivity.  At  scales  larger  than 
a  meter,  airfoil  data  are  unreliable  owing  to  increasing  vehicle  motion  and  uncertain 
and  varying  probe  response. 


ocean  velocity  profiles.  The  ACM  results  are  important  in  that  they  can  resolve 
shear  variance  produced  by  internal  waves  and  turbulence.  To  interpret  the  relative 
ACM  measurements,  we  model  the  motion  of  the  vehicle  as  it  responds  to  oceanic 
shear.  The  tail  array  of  drag  brushes  and  turning  blades  causes  MSP  to  respond 
more  severely  than  other,  similar  profilers  such  as  TOPS  and  HRP.  Modification  of 
previous  response  models  of  H  ,ves  et  al.  (1984)  and  Schmitt  et  al.  (1988)  to  include  a 
proper  specification  of  the  large  tail  force  was  a  major  advance  in  the  MSP  analysis. 
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In  Section  II,  we  describe  the  design  and  operation  of  MSP,  discuss  its  behavior 
as  a  platform,  and  anticipate  eifects  on  the  ACM  measurements.  In  Section  III 
we  examine  the  sensors  and  their  data,  emphasizing  the  transformation  of  ECM 
and  ACM  measurements  into  earth  coordinates.  In  Section  IV,  after  considering 
vehicle  kinematics  and  dynamics,  we  present  our  technique  for  converting  the  ACM 
measurements  into  ocean  velocity  profiles — correction  for  tilt  and  for  attenuation  of 
signals  longer  than  MSP  due  to  its  gross  motion  as  a  point  mass.  We  explain  how  we 
use  analytic  transfer  functions  to  find  optimal  values  for  the  model  parameters.  In 
Section  V,  we  scrutinize  ocean  velocity  profiles  and  shear  spectra  from  the  ACM  and 
ECM  (and  airfoil  probes)  and  discuss  the  strengths  and  weaknesses  of  the  results.  In 
the  final  section,  we  summarize  the  important  aspects  of  MSP  behavior  that  affect  its 
ability  to  resolve  oceanic  shear  variance  over  a  wide  band  of  scales.  Five  appendices 
supplement  this  report.  Appendix  A  provides  greater  detail  about,  the  sensors  and 
processing  of  their  data.  Appendix  B  deals  with  the  treatment  of  airfoil  probe  data 
and  compares  airfoil  and  ACM  spectra.  Appendix  C  discusses  the  vertical,  horizontal, 
and  till  dynamics  of  the  vehicle.  Appendix  D  details  the  numerical  scheme  employed 
for  the  response  model  integration.  Appendix  E  derives  the  model  transfer  functions 
a  -d  uses  them  to  demonstrate  that  tilt  dominates  the  accelerometer  data. 
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II.  VEHICLE 

A.  Vehicle  Design  and  Operations 

Much  of  the  design  of  MSP  is  aimed  at  providing  a  suitable  platform  for  its 
sensors.  As  an  untethered  and  freely  falling  vehicle,  it  largely  avoids  high-frequency 
vibrations  and  oscillations  that  could  confound  measurement  of  turbulent  velocity 
fluctuations.  It  descends  slowly  and  steadily  for  adequate  resolution  of  temperature 
microstructure,  and  it  spins  about  its  axis  for  reliable  measurement  of  ambient  electric 
fields. 

Most  of  the  length,  mass,  and  volume  of  MSP  comes  from  its  main  pressure 
case,  which  houses  the  batteries,  microcomputer,  memory,  electronics,  and  vehicle- 
orientation  sensors — accelerometers  to  monitor  vehicle  tilt  and  vibration  and  a  magne¬ 
tometer  to  monitor  rotational  orientation  (Figure  1).  A  gap  of  a  few  inches  separates 
the  white  outer  skin  from  the  pressure  case.  Syntactic  foam  near  the  tail  fills  much  of 
this  space,  and  seawater  floods  the  rest  upon  deployment.  The  outer  skin  streamlines 
the  instrument  and  forms  a  nonconducting  cylinder  between  the  electrodes  of  the 
ECM.  To  descend  while  acquiring  data  and  then  to  ascend  for  recovery,  MSP  uses 
releasable  ballast  weights,  without  which  its  displacement  renders  it  positively  buoy¬ 
ant  (Table  2).  Attached  to  the  tail  are  drag  brushes  to  slow  the  fall  and  adjustable 
turning  blades  to  set  the  rotation  rate. 

The  leading  end  (nose)  of  MSP  is  home  to  several  sensors  and  the  two  ballast 
tubes  (Figure  2).  A  guard  cage  protects  the  sensors.  The  ACM  is  the  first  to  sample 
at  a  given  depth,  as  its  transducers  are  positioned  below  the  end  of  the  cage  to 
eliminate  spurious  reflections.  Next  comes  a  cluster  of  fine-  and  microstructure  sensors 
(Table  3),  positioned  to  sample  fluid  unperturbed  by  the  ACM  or  cage.  Above  the 
sensors  a  spherical  plastic  cap  rounds  the  blunt  face  of  the  instrument. 

Each  deployment  begins  with  the  operator  recharging  the  batteries  and  setting 
several  drop  parameters.  External  cables  are  disconnected,  and  MSP  is  put  into 
the  water.  Temporary  floats  (beach  balls)  hold  the  instrument  at  the  surface  while 
the  ship  moves  away  to  avoid  contamination  of  the  ECM  measurements.  After  a 
few  minutes  MSP  releases  the  beach  balls  and  starts  to  fall.  During  descent,  data 
are  recorded  in  20  Mbytes  of  solid  state  memory;  this  storage  is  volatile  in  that 
data  are  lost  if  the  batteries  lose  their  charge  (typically  aft  -  S  hours).  When  the 
microcomputer  determines  that  either  the  pressure  or  time  limit  has  been  reached, 
the  ballast  is  released  to  end  the  drop.  Corrosible  links  and  fracture  pins  back  up 
the  primary  mechanisms,  and  in  shallow  waters  an  acoustic  detector  aborts  the  drop 
before  striking  the  bottom.  While  ascending,  MSP  transmits  an  acoustic  telemetry 
signal  containing  its  slant  range  and  depth  (similar  to  the  procedure  described  by 
Sanford  et  al.,  1978)  so  that  the  ship  can  approach  the  rising  instrument.  At  the 
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Figure  1.  Diagram  of  Multi-Scale  Profiler  and  photograph  of  post-drop  recovery. 
Length  scale  is  for  the  diagram  at  left.  Blades  and  brushes  are  drawn  for  the  large- 
tail  configuration;  dashed  lines  indicate  width  of  the  small  tail.  Positions  of  centers 
of  buoyancy  and  mass,  Cy  and  C 'a .  are  approximate,  owing  to  uncertainties  and  to 
variations  in  ballast  mass  and  density  of  trapped  water.  The  pressure  sensor,  hidden 
among  the  turning  blades,  is  attached  atop  the  main  pressure  case.  Photo  at  right  shows 
recovery  after  a  Florida  Straits  drop  with  the  small  tail. 
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Table  2.  Vehicle  parameters.  Rightmost  column  indicates  notation  used  in  the  main 
text.  Moments  are  for  transverse  rotations  relative  to  the  ECM  (between  CG  and  CB ). 
Values  in  parentheses  are  for  the  small-tail  array,  if  different  than  for  the  large  array. 


Total  length 

4.3 

m 

Radius,  main  tube 

0.21 

m 

(r/v) 

Radius,  tail  array 

0.57  (0.42) 

m 

(rr) 

Pressure  tube,  length 

2.6 

m 

Pressure  tube,  radius 

0.16 

m 

Vehicle  displacement 

0.27 

m3 

(r.) 

Vehicle  mass,  with  ballast 

280 

kg 

(M.) 

Ballast  mass 

10 

kg 

Trapped  water  volume 

0.10 

m3 

Total  mass,  including  trapped  water 

382 

kg 

(AO 

Center  of  mass,  from  ACM 

1.90 

m 

(Co) 

Center  of  buoyancy,  from  ACM 

2.05 

m 

(CB) 

Moment  of  inertia,  effective 

830  (730) 

kg  m2 

(Ie) 

Damping  coefficient,  effective 

980  (610) 

kg  m2  s  1 

(Re) 

Righting  moment,  effective 

580  (550) 

N  m 

(Re) 

Maximum  pressure 

12.0 

MPa 

Effective  added  mass  displacement 

0.40 

m3 

(re) 

Water  density,  nominal 

1025 

kg  m-3 

(p) 

surface,  the  radio  and  flashing  light  at  the  tail  (and  the  bright  orange  drag  brushes) 
help  visual  sighting.  A  small  boat  (e.g.,  Zodiac)  is  launched  to  gather  MSP  and  drag 
it  safely  to  its  pickup  point  beside  the  ship  (Figure  1).  Once  aboard,  the  instrument 
is  hooked  up  to  download  its  data,  and  preparations  are  begun  for  the  next  drop. 
Total  turnaround  time  for  a  10-MPa  (1000-m)  drop  is  around  3  hours. 

B.  Vehicle  Motion 

Table  1  lists  nominal  fall  and  rotation  rates  for  the  five  major  MSP  cruises. 
Although  faster  rotation  would  enhance  ECM  performance,  the  associated  dynamics 
could  complicate  the  ACM  and  airfoil  measurements.  Intercruise  variations  in  fall 
and  spin  rate  generally  reflect  alterations  in  the  disposable  ballasl  and  blade  angles. 
We  switched  to  the  small  tail  configuration  in  1990.  which  weakened  lateral  forcing 
to  improve  the  ACM  measurements  but  also  reduced  vertical  drag  thereby  increasing 
the  fall  rate. 
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Figure  2.  Photograph  of  MSP  nose  .section  with  sensors.  From  top  lc>  bottom,  sensors 
in  the  central  cluster  arc1:  Noil  Brown  conductivity,  airfoil  probe  1.  FP07  temperature 
( linesi  met  lire),  airfoil  probe  2.  and  FP07  temporal  tin'  (gradient ).  For  clarity,  cables  and 
pump  hose  are  disconnected  Irom  the  Sea-Bird  sensors.  Second  ballast  tube  is  obscured 
bill  is  diametrically  opposed  to  the  first.  Coiard  rings  (center  and  right)  are  separated 
by  (1.2b  in. 


Mu’  MSP  falN  slowly.  l).2">  D.-'ib  in  s~‘.  to  minimize  attenuation  by  the  glass 
bead  thermistor  ol  dissipal ion-scale  temperature  lliiet nations.  We  attain  a  faille 
uniform  descent  by  ballasting  for  around  d.b-kg  negative  buoyancy,  a  \  hie  that  i- 
large  compared  with  variations  in  buoyancy  that  occur  during  a  drop  as  increasing 
water  density  competes  with  shrinking  ol  MSP  via  cooling  and  compression.  Owing 
to  the  broadened  axial  cross  section  due  to  the  ring  of  polyethylene  brushes  at  the 
tail,  the  vertical  quadratic  drag  balances  the  negative  buoyancy  at  acceptably  low 
bill  rales.  |  |u.  brushes  also  homogenize  the  wake  and  st i lie  vibration-  that  might 
arise  t  lirough  eddy  shedding,  kali  rate  profiles  (  h  igure  3)  show  sonic  consistency  wit  h 
the  simple  buoyancy-drag  balance:  in  the  lirst  II)  20  rn.  MSP  accelerates  toward  a 
terminal  velocity;  as  its  buoyancy  increases  in  tin'  sharpest  part  of  the  pyenocline. 
MSP  slows  rapidly  (wit  h  an  ext  reme  (IW'fJl  —  ().()()  I  in  s“*  in  drop  02 1  I  at  ]  MPa): 
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Table  3.  Sensor  specifications. 


Sensor 

Sample 

Rate 

(Hz) 

lypeof  No. 
Digiti/a-  of 
tion*  Axes 

Range 

Resolution 

Noise, 
nns  b 

Magnetometer 

25 

Fib  3 

±60000  nT 

30 

10 

degrees,  rotation  or  lilt:r 

0.05 

0.07 

Accelerometers 

25 

16b  2 

±16.4  in  s ~2 

5.0  x  10“4 

2.6  x  lO"4 

degrees  tilt: 

0.003 

0.0015 

ECM  d 

25 

ftc  2 

±150  rn  s"1 

1.0  x  10~a 

2.0  x  10-3 

ACM 

62.5 

16b  2 

±0.83  m  s“: 

2.5  x  I0“5 

4.2  x  10~4 

Airfoil  Probes  f 

125 

16b  2 

±49  s"1 

1.5  x  10“3 

(App  B) 

Pressure 

25 

ftc 

0  21  MPa 

3.9  x  10-“ 

1.8  x  10~4 

Temp.,  Sea-Bird 

25 

ftc 

1  31  °C 

1.6  x  10~4 

2.3  x  10"4 

Cond.,  Sea-Bird 

25 

ftc 

1  -  6  S in-1 

3.2  x  10-5 

2.8  x  10'5 

Temp.,  FP07 

125 

16b 

0-35  °C 

8.0  x  10-“ 

4.0  x  10-4 

Cond.,  Neil  Brow-n 

125 

16b 

3  -  6.2  S  m~J 

5.0  x  10“5 

2.5  x  10~5 

Temp.  Grad.,  FP07 

f  250 

16b 

±18  °C  in"1 

5.4  x  10~4 

— 

°12b  =  12-bit  A/D.  ±10  V  input;  16b  =  16-bit.  A/D,  ±8.192  V  input;  ftc  =  frequency- 
time  counter,  two  16-bit  words  output  (frequency  range  of  input  signal  depends  on  sensor: 
ECM.  5-10  kHz;  pressure,  34-40  kHz;  Sea-Bird,  6-11  kHz). 

6Noise  in  raw  data,  for  a  single  component  (axis),  computed  by  integrating  the  spectral 
noise  floor  out  to  the  Nyquist  frequency. 

cComputed  with  nominal  values  of  30,000-50,000  nT  for  components  of  the  earth’s  magnetic 
field. 

dECM  velocities  vary  inversely  with  Bz .  the  vertical  component  of  the  geomagnetic  field; 
listed  values  are  for  Bz  =  40, 000  nT. 

cThe  two  axes  of  the  airfoil  probes  are  usually  aligned  in  parallel  for  redundancy;  range 
and  resolution,  valid  for  1-20  cpm  and  0.30  m  s-1,  increase  with  higher  w-avenumber  or 
slower  fall  rate;  dissipation  rates,  e,  computed  from  airfoil  probe  spectra  have  a  noise  level 
of  (3-10)xl0~n  W  kg'1. 

Range  and  resolution,  valid  for  1-7  cpm,  increase  with  higher  wavenumber  (some  temper¬ 
ature  variation,  also);  dissipation  rates,  \,  computed  from  gradient  spectra,  have  a  noise 
level  of  ( 1  ()— 1 2  -  10_H)  °C2  s'1. 
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Figure  3,  Profiles  of  fall  rate;-  (heavy  lines)  ami  in  silu  density  (light  lines)  for  dro[>s 
during  three  MSP  cruises.  Subsequent  are  offset  by  0.1  in  s*1.  rr,tp  by  K  kg  i»-\ 

Below  the  top  of  the  pycnodine,  H  in,,,  =  (0.31  ±0.02)  in  s-1  for  0241.  (0.275  ±  0.02) 
for  0211.  and  (0.25  ±  0.01 )  for  0000. 


however,  the  fall  rate  shows  no  clear  tendency  thereafter.  Large-scale  fluctuations  of 
0.01  0.03  rn  s-1  are  likely  due  to  vertical  motions  of  the  water,  such  as  internal  waves, 
carrying  MSP  along  with  little  effect  on  the  relative  (low  (l)esaubies  and  Gregg,  1978). 
DespUe  uncertainties  in  net  buoyancy  and  relative  flow,  we  estimate  a  vertical  drag 
coefficient  of  1.0±0.15  based  on  comparison  of  tall  and  rise  rales  just  before  and  after 
ballast  release. 

Currents  passing  across  MSP  induce  quadratic  drag  and  linear  pressure  forces 
which  accelerate  the  vehicle  horizontally.  Dominating  the  response  are  the  pressure 
forces,  which  occur  at  the  nose  and  tail,  where  the  vehicle  radius  widens.  The  nature 
of  MSP's  response  to  velocity  fluctuations  depends  on  their  vertical  scale.  Owing 
to  the  nose  force  and,  to  a  lesser  extent,  lateral  drag  on  the  main  tube,  MSP  fully 
tracks  the  large-scale  flow,  although  its  inertia  causes  it  to  lag  behind  somewhat.  The 
tail  force  quickens  MSP’s  response  such  that  it  can  partial!}'  follow  fluctuations  at 
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near  !«•  «.<  al«->  Ithereby  romplii  at  mg  interpretation  of  the  ACM  measurements), 
i  In-  <1im iiifpiushes  MSI*  from  less  responsive  profilers  such  as  TOPS  (Hayes  et  al., 
19Ml  and  HUP  (Schmitt  et  al..  1988).  which  have  no  comparable  tail  structure. 

In  its  tilt  behavior.  MSP  acts  like  a  furcitl,  dampt'd  pendulum.  Separation 
bet  "<vn  the  centers  of  buoyancy  an«l  mass  (including  the  trapped  water;  estimated 
in  Iable  2)  provides  most  of  the  righting  moment,  and  drag  along  the  main  tube 
ri'sists  oscillations.  As  MSP  falls,  further  contributions  arise  from  forces  at  the  tail: 
vertical  drag  adds  to  the  righting  moment,  and  oscillation-induced  pressure  forces 
produce  strong  damping.  The  nose  and  tail  forc«*s  from  the  lateral  flow  generate  a 
forced  response,  which  ovei whelms  free  oscillation.  Prom  incidents  in  which  MSP 
was  impulsively  forced  (by  shark  hits?),  we  observe  a  natural  period  at  8-11  s  being 
dampen!  out  within  two  oscillations,  consistent  with  the  moments  estimated  in  Table  2. 

The  effect  of  MSP's  horizontal  response  on  the  ACM  depends  on  the  scale  of 
the  fluctuation  being  measured.  To  focus  on  the  gross  vehicle  motion,  we  disregard 
rotation  and  tilt.  At  depth  c  (positive-downward),  the  ACM  measures  the  relative 
motion 

l'»ciii(~)  =  T’w«i(*!  T’m»j>(“) 

between  the  water  and  vehicle.  Casting  the  time  lag  of  MSP  as  a  distance  lag,  D , 
from  a  velocity  fluctuation  to  the  ACM  position,  we  express  the  vehicle  motion  in 
response  to  a  given  wavelength.  A,  as 

'W(')  =  7(^)»’w«i ('-/))  , 

where  -)  is  the  response  factor;  7  =  1  for  full  response  at  large  scales  and  0  for  no 
response  at  small  scales  (since  inertia  prevents  MSP  from  reacting  to  rapid  fluctua¬ 
tions).  The  ACM  measurement  becomes 

t Wm  (  -  )  T-’wa[  (  C  )  7  L’wat  (  -  D)  » 

which  reduces  to  iv-m(')  =  tv«i ( - )  —  t’w« (:  —  D)  for  large  scales  (such  as  constant 
shear  flows)  and  to  i>*cni(-)  =  t'«ii(-)  for  small  scales.  The  strong  tail  force  causes 
MSP  to  respond  (7  >  0)  at  scales  near  its  own  length  and  reinforces  the  nose  force  to 
establish  a  lag  of  about  one  instrument  length.  In  contrast,  TOPS  and  HRP  respond 
little  (7  ss  0)  at  scales  smaller  than  10  m  and  have  lags  closer  to  two  instrument 
lengths.  Figure  4  illustrates  important  tendencies  in  the  ACM  measurements  for 
various  scales  of  water  motion,  uwat(r)  ~  sin(2n/\)z .  At  A  >■  D  the  measurement 
approaches  D  </ewat/ dz,  the  constant  shear  limit.  The  ACM  moves  against 

currents  of  A  —  2D  and  with  those  of  A  =  D,  measuring  velocities  respectively 
amplified  by  (1  +  7)  and  diminished  by  (1  —7).  These  effects  complicate  our  measured 
ACM  spectra,  yet  hardly  perturb  spectra  from  the  less  responsive  TOPS  and  HRP. 
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10D  (y=  1) 


Figure  4.  Effects  of  MSP  point-mass  response  on  A('M  measurements  due  to  water 
motions  of  various  wavelengths,  X.  Heavy  curves  are  water  velocity.  t;w»t  ~  sin  A”,  where 
k  —  ‘Jtt / X  and  c  is  depth,  positive-downward.  Light  curves  are  ACM  measured  velocity, 
IV, ,i  =  rw*t  -  For  the  vehicle  velocity.  enu|) ,  the  response  factor,  7,  lessens  with 

scale  (see  text );  for  this  illustration,  the  response  lag  distance  of  D  -  5  m  slightly  exceeds 
MSP's  length.  At  large  scales,  »Vui  — 1 ■  D  dvva,i/dz.  At  A  =  2 D,  tvni  =  (1  +  7)ew*i'.  at 
A  —  V .  i-Vm  =  (1  -  7)ew*i  (inflated  7’s  are  used  here  to  clearly  illustrate  responses). 
When  7  =  0  at  small  scales,  the  ACM  measures  the  water  directly,  tv™  =  rw>t.  Not 
considered  here  are  contributions  to  ACM  data  by  vehicle  oscillation. 


Also,  tilt  orientation  and  oscillation  affect  ACM  measurements  in  a  manner  that 
intensifies  our  spectral  complications. 

Before  considering  further  the  vehicle  dynamics  and  associated  ACM  effects,  we 
investigate  the  array  of  sensors  and  the  measurements  that  they  obtain. 
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III.  MEASUREMENTS 

Table  3  contains  specifications  for  the  MSP  sensors.  Range  and  resolution  de¬ 
pend  on  tin*  digitization  method,  as  signals  are  passed  to  either  an  A/D  converter 
or  a  frequency-time  counter.  Conversion  from  voltage  or  frequency  to  scientific  units 
employs  calibrations  based  on  the  static  response  of  the  sensor  and  its  associated 
electronics.  Noise  levels  are  estimated  from  spectra  of  data  over  quiet  regions  of  the 
water  column:  the  spectral  noise  floor  is  integrated  out  to  the  Nyquist  frequency 
(Appendix  A).  Vertical  resolution  is  limited  by  noise  and  also  by  the  dynamic  re¬ 
sponse  associates!  with  sensor  size  or  response  time.  For  those  devices  that  measure 
along  multiple  directions.  Table  3  pertains  to  a  single  component;  an  individual  data 
channel  corresponds  to  each  sensitive  axis.  Axes  are  defined  relative  to  an  orthogo 
nal  reference  frame,  ijk,  fixed  to  MSP.  In  Appendix  A,  we  discuss  the  mapping  of 
data  from  sensor  frames  onto  tin*  common  ijk  frame  and  subsequently  to  the  eart  h's 
xyz  frame. 

Fur  each  drop,  smoothed  pressure  data  are  used  to  build  an  index  file  that 
indicates,  for  a  specified  pressure  and  data  channel,  the  elapsed  time  at  which  a 
measurement  was  taken,  (juided  by  the  index,  raw  data  are  low  pass  filtered  and  then 
subsampied  onto  evenly  spaced  pressure  grids  (Table  4).  Reference  grid  pressures  are 
also  used  to  estimate  fall  rates  and  to  compute  hydrographic  quantities. 


A.  CTD 

Two  pairs  of  conductivity-temperature  sensors  provide  redundant  measurements 
ami  enable  intercomparison  of  performance  and  response.  The  Sea-Bird  pair  has 
better  stability,  while  the  faster  response  of  the  other  pair,  a  Neil  Brown  conductivity 
cell  and  a  Thermometries  Fast  ip  FP07  thermistor,  provides  finer  vertical  resolution. 
As  of  1W9U.  ducting  between  the  Sea-Bird  sensors  ensures  that  fluid  pumped  through 
the  conductivity  cell  is  first  sampled  by  the  thermistor.  We  calibrate  the  Neil  Brown 
conductivity  cell  and  Thermometries  Fast  ip  FP07  thermistor  against  the  standard  of 
the  Sea  Bird  pair,  using  either  in  situ  drop  data  or  measurements  in  a  controlled  bath 
in  our  laboratory. 

Results  from  the  CTD  data  include  profiles  of  temperature,  salinity,  density,  and 
buoyancy  frequency  and  vertical  wavenumber  spectra  of  temperature  and  displace¬ 
ment  (or  strain). 

B.  Microstructure 

Turbulent  fluctuations  in  temperature  are  sensed  by  a  second  FP07  thermistor. 
Its  signals  are  passed  through  a  differentiating  circuit  to  emphasize  the  smaller  scales 
and  a  60- Hz  four- pole  Chebyschev  filter  to  suppress  high-frequency  noise.  As  we  form 
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Table  4.  Data  processing  specifications  for  filtering  and  subsampling.  Grid  spacing  is 
for  standard  processing;  in  seawater,  a  pressure  interval  of  10  kPa  =  0.01  MPa  closely 
corresponds  to  a  depth  interval  of  1  m.  For  the  filter,  when  “BW"  is  appended,  the 
full  length  of  the  low-pass  Bartlett  window  is  indicated;  for  the  microstructure  data, 
parentheses  enclose  the  sizes  of  spectral  processing  windows  and  segments  (w,s)  (sec1 
Appendix  B). 


Sensors 

Grid 

Filter 

Comments 

Pressure® 

0.96  s 

CTD 

1.92-s  B\Y 

Interpolated  onto  10-kPa  grid 

Sea  Bird 

1  kPa 

2-kPa  BYY 

Ducted 

5  kPa 

20-kPa  B\Y 

Not  ducted 

Neil  Brown  and  FPUT 

1  kPa 

1-kPa  BW 

Microstructuir  Sped 

ra 

Temperature  ( \ ) 

20  kPa 

(24, G)  kPa 

Robust 

5  kPa 

(5,5)  kPa 

Periodogram 

Airfoil  probes  (< ) 

20  kPa 

(24,6)  kPa 

Robust 

5  kPa 

(5,5)  kPa 

Periodogram 

KCM  '■ 

5  kPa 

\'elocity 

10-kPa  BW 

ACM' 

1  kPa 

2-kPa  BW 

“Smoothed  pressure  data  are  subsampled  onto  a  temporal  grid,  from  which  clasped  drop 
times  are  interpolated  every  10  kPa  to  form  the  pressure  reference  index. 

^'Afler  being  oriented  via  magnetometer  and  accelerometer  data,  drift-filtered  ECM  mea¬ 
surements  are  converted  to  velocity  using  the  geomagnetic  field,  B  -  ll^y  +  Hz  z. 

f  After  being  oriented  via  magnetometer  data,  ACM  relative  measurements  are  converted  to 
ocean  velocity  by  using  accelerometer  data  and  the  point -mass  response  model  to  remove 
platform  motion. 


power  spectra  from  these  data,  we  account  for  the  static  calibration  and  dynamic 
response  of  the  probes  (Gregg  and  Meagher,  1980)  along  with  the  effects  of  the 
electronics.  The  spectra  are  used  for  comparing  microstructure  with  finestructure 
and  to  estimate  \. 

The  pair  of  airfoil  probes  detects  fluctuations  in  horizontal  velocity  (Osborn  and 
Crawford,  I960).  For  redundancy,  we  align  both  probes  in  the  same  direction.  As 
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the  water  pushes  oil  the  probe  tip,  the  deflection  of  an  enclosed  piezobimorph  beam 
generates  a  voltage,  which  is  then  modified  by  a  differentiator  and  a  30-Hz  four-pole 
Chebvschev  filter.  The  probe  tip  smooths  out  fluctuations  of  scale  smaller  than  2  cm 
(Ninnis,  1984).  We  divide  raw  data  spectra  by  the  combined  transfer  function  of  the 
probe  and  electronics  to  produce  microscale  shear  spectra,  which  we  then  integrate 
to  estimate  i  (Appendix  B).  By  patching  ensemble- averaged  spectra  onto  those  from 
the  ACM,  we  broaden  the  bandw-idth  to  100  cpm  and  cross  check  spectral  levels  in 
the  overlapping  band.  Airfoil  spectra  are  unreliable  at  wavenumbers  below  1  cpm. 
owing  to  vehicle  motion  and  a  poorly  understood,  probe-dependent  response  involving 
temperature  sensitivity. 

C.  Vehicle  Orientation 

A  Develco  fluxgate  magnetometer  monitors  the  rotational  position  and  motion 
of  the  vehicle.  The  horizontal  data  allow  transformation  of  ECM,  ACM,  and  ac¬ 
celerometer  measurements  from  the  instrument  frame  to  geomagnetic  coordinates. 
The  data  reflect  not  only  the  geomagnetic  field,  B,  but  also  the  magnetic  signature 
of  the  instrument  and  effects  of  uncertain  sensor  alignment,  bias,  and  sensitivity. 
With  the  value  of  B  at  a  drop  location  approximated  by  a  spherical  harmonic  model, 
the  nongeomagnetic  contributions  can  be  estimated  and  removed  from  the  data  (Ap¬ 
pendix  A). 

Sundstrand  Q-flex  linear  servo- accelerometers  are  aligned  along  three-  orthogonal 
axes,  hut  only  the  two  horizontal  components  are  used.  Their  primary  function  is 
to  monitor  vehicle  tilt,  which  is  the  dominant  signal  (Appendix  E).  Spectra  hit  a 
flat  noise  floor  at  1  to  2  cpm;  the  lack  of  spectral  spikes  above  this  low-level  floor 
indicates  that  there  are  no  high-frequency  vibrations  strong  enough  to  contaminate 
the  airfoil  data.  Although  accelerometer  data  can  augment  the  magnetometer  data 
for  more  precise  determination  of  rotational  orientation,  the  effect  is  inconsequential 
for  small  tilt  angles.  On  the  other  hand,  the  variation  in  tilt  demonstrates  oscillation 
and  wobbling  of  the  platform,  motions  that  significantly  affect  the  ACM. 


D.  Velocity 

The  ECM  detects  voltages  induced  by  ocean  currents,  a  technique  nearly  insen¬ 
sitive  to  horizontal  motion  of  the  vehicle.  Unfortunately,  thermal  gradients  between 
the  electrodes  lead  to  large  drifting  offsets  in  the  data.  We  isolate  the  oceanic  signal 
by  using  a  low-frequency  spline  to  estimate  the  drift.  However,  the  results  are  suspect, 
(especially  for  subrotation  scales)  because  the  drift  is  often  comparable  to  the  signal 
amplitude  over  a  rotation  (Figure  5).  With  our  two  pairs  of  electrodes,  we  could  gen¬ 
erate  velocities  independently  from  each  pair  as  done  by  Sanford  et  al.  (1978),  and 
then  compare  or  average  the  results;  vertical  resolution  is  on  the  order  of  the  rotation 
scale  (A'~,\  where  A'lol  =  Qrot/ir  is  the  rotation  wavenumber  in  cycles  per  meter). 
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figure  5.  Data  from  i  and  j  ECM  axes.  In  tin*  left  panel.  ECM  measurements  arc' 
ctd.M'i  (mc'iusun'd  c ,  and  c arc,  respectively,  1100  and  ISO  /<  V  at  (i  MPa)  and  wrapped 
around  (for  tin*  high-drifting  t  axis).  Differences  between  measurements  (heavy)  and 
estimated  ollsc'ts  (light)  arc  the  oceanic  signals,  which  are  shown  in  the  right  panel  (c, 
is  hea\  y;  < )  is  light ). 


Instead,  we  compute  the  velocity  at  each  grid  point  In  interpreting  the  filtered  data 
as  orthogonal  components  of  the  eleciric  field  vector  (Appendix  A).  This  method 
potentially  yields  liner  resolution,  hut  the  useful  limit  remains  near  the  rotation  scale 
owing  to  uncertainty  in  the  removed  drift.  Since  ECM  sensitivity  is  proportional  to 
/f/,  the  vertical  component  of  B,  performance  degrades  near  the  magnetic  equator 
with  diminishing  signal  to  noise  ratio.  In  spite  of  their  restrictions.  ECM  results  arc- 
critical  in  gauging  the  low- wavenumber  performance  of  the-  vehicle  response-  model. 

The  Neil  Brown  Instrument  Systems  ACM  measures  the  relative  flow  along  the 
0.2-m  path  between  transducers.  Two  transducer  pairs  define  orthogonal  horizontal 
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axes.  T  heir  low  noise  level  (Table  3)  allows  the  ACM  to  resolve  motions  down  to 
vertical  scales  less  than  a  meter.  Because  of  MSP’s  motion,  the  ACM  detects  a 
velocity  different  than  that  of  the  water  (as  reflected  in  Figure  4).  As  such,  our 
ability  to  generate  ocean  velocities — particularly  at  scales  greater  than  a  few  meters— 
is  constrained  by  the  success  of  the  vehicle  response  model.  At  high  wavenumbers, 
measurements  are  attenuated  for  motions  whose  horizontal  scales  approach  the  0.2- 
m  path  length.  Also,  small-scale  data  are  intermittently  contaminated  by  narrow- 
banded  oscillations  of  drifting  frequency,  which  appear  as  a  wandering  spike  in  ACM 
spectra  (Appendix  A). 

For  both  ECM  and  ACM  processing,  measurements  must  be  transformed  into 
geomagnetic  coordinates.  In  preparation,  velocity  and  navigation  data  must  be  cor¬ 
rected  for  bias  and  axis  alignment  (Appendix  A).  Once  transformed,  ECM  data  can 
be  directly  converted  into  ocean  velocity  profiles,  but  ACM  data  cannot  since  they 
remain  affected  by  platform  motion  and  tilt. 

In  Figure  6,  we  compare  ECM  and  ACM  profiles  at  this  stage  of  processing  for 
a  portion  of  PATCIIEX  drop  0(J9(i.  We  will  follow  the  development  of  this  ACM 
record  into  an  ocean  velocity  profile.  Table  5  presents  nomenclature  relevant  to  these 
processed  data  and  to  subsequent  processing  in  the  vehicle  response  model.  The 
oceanic  velocity  estimated  from  the  ECM  is  denoted  vem,  while  that  resulting  from 
application  of  the  response  model  to  the  ACM  measurements,  jv,,,,  is  denoted 

For  the  profiles  in  Figure  6,  maximum  relative  flows  at  the  ACM  and  vehicle  tilts 
are  about  0.03  m  s-1  and  0.3°.  In  data  from  other  cruises,  relative  speeds  reach  over 
0.05  m  s_l  at  sharp  shear  features  and  up  to  0.10-0.20  m  s_1  within  zones  of  high 
large-scale  shear  (such  as  above  the  Equatorial  Undercurrent).  Tilt  angles  seldom 
exceed  1°,  except  when  influenced  by  surface  waves  or  strong  shear. 

Figure  7  displays  power  spectra  of  ECM,  ACM,  and  accelerometer  measure¬ 
ments  from  28  PATCIIEX  drops  (including  0090).  These  are  ensemble  averages  of 
total  spectra  (east  and  north  autospectra  combined)  from  1-MPa  (100-m)  segments 
spanning  4-8  MPa  (further  details  in  Appendix  A).  Cregg  and  Sanford  (1988)  report 
shear  levels  during  PATCHEX  comparable  to  those  of  the  Garrett  and  Munk  (GM70) 
internal  wave  model  (Garrett  and  Munk,  1975;  see  also  Gregg  and  Kunze,  1991).  In¬ 
deed,  the  ECM  spectrum  is  dose  to  CM7(i  at  wavenumbers  lower  than  Ar„.;  spectral 
features  near  Ar<„  and  2  Ari(t  (and  the  divergence  from  GM7C  at  higher  wavenumbers) 
reflect  shortcomings  of  the  ECM  processing  related  mainly  to  the  large  signal  drift. 
The  ACM  spectrum  fluctuates  about  the  ECM  (and  CM7G)  spectrum,  displaying  fea¬ 
tures  suggested  in  Figure  4:  it  has  excess  variance  at  0.1  >  pm  ( A  2  D)  and  a  deficit 
near  0.2  cptn  (A  ss  D)  which  we  refer  to  as  the  u0.2-cpm  notch.”  When  corrected 
for  the  variance  removed  when  data  were  smoothed  with  the  2-kPa  Bartlett,  window, 
the  ACM  spectrum  is  well  above  the  noise  floor.  Past  3.5  cpm,  aliasing  contaminates 
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Figure  6.  PATCH  EX  drop  009G,  profiles  of  measurements  after  transformation  to 
east  (left)  and  north  (right)  components.  ECM  velocities  (light  lines)  are  demeaned 
for  display;  ACM  velocities  (heavy  lines)  are  not.  Tilt  angles,  8,  are  computed  from 
accelerometer  data  as  8  —  (anieas/y),  in  radians.  Note  that  1  MPa-  100  m  relates 
pressure  to  depth,  approximately. 


the  spectrum  (Appendix  A).  The  broad  peak  of  the  accelerometer  spectrum  near 
0.12  cpm  reflects  the  forced  tilt  response  as  velocity  fluctuations  push  the  nose  and 
tail  in  opposite  directions. 

To  generate  an  oceanic  spectrum  from  the  ACM  measurements,  we  must  correct 
them  for  tilt  and  for  MSP’s  gross  motion  as  a  point  mass.  We  seek  to  convert  the 
distorted  measurement  spectrum  into  the  true  oceanic  spectrum.  We  require  that 
the  final  spectrum  be  reasonably  smooth  and  that  it  be  close  to  ECM  levels  at  low 
wavenumbers.  With  the  PATCHEX  dataset  as  a  lest  case,  the  GM7b  model  spectrum 
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can  also  provide  some  guidance.  To  reduce  the  dynamic  range  and  emphasize  spectral 
features,  subsequent  displays  are  in  the  form  of  shear  spectra,  computed  as  $a{k)  = 
(2tt k)2$v(k),  with  wavenumber  k  in  cycles  per  meter. 

With  ACM  and  accelerometer  measurements  in  geomagnetic  coordinates,  we 
assume  that  we  have  accounted  for  all  rotational  effects.  We  have  evaluated  spin- 
related  dynamics,  such  as  gyroscopic  and  Magnus  effects,  and  found  them  negligible 
relative  to  those  included  in  our  response  model.  Therefore,  we  can  consider  motions 
and  forces  in  the  xz-plane  independently  of  those  in  the  1/2-plane. 


Table  5.  Measured  or  derived  variables  in  velocity  processing.  Functions  of  depth,  2,  are 
attributes  of  the  water,  independent  of  MSP.  Those  describing  or  involving  vehicle  motion  arc 
defined  as  functions  of  t .  the  elapsed  time  of  the  drop  when  the  ACM  is  at  depth  2(0.  Fourier 
transforms,  e.g.,  fiwot,  are  defined  in  Appendix  E. 


Variable 

Description 

Spectrum 

Transfer  Function 

[m  s~‘] 

Water  velocity 

*  V 

1 

*w(0 

[in  s-1] 

Vehicle  velocity  (as  point  mass) 

$«*»!• 

A{k)  = 

M'nap(0 

[in  s_1] 

Absolute  fall  rate 

0(t) 

[radians] 

Vehicle  tilt 

B(k)  = 

0  /  t)wal 

s{t) 

[m] 

Depth  at  ACM 

t’rnl ( 0 

[ill  s  *] 

Relative  horiz.  vel.  (at  ACM) 

<f,re! 

V 

Cr(k)  = 

t’acn.  (  0 

[m  s_I] 

ACM  measured  velocity 

V 

Cm(k)  = 

T’ttcjii  /  f’wat 

^meas(  ^  ) 

[m  s-2] 

Measured  acceleration 

(fttneas 

D(k)  = 

^meas  /  ^Vat 

VN(t) 

[ms  *] 

Relative  velocity  at  nose 

VT(t) 

[in  s"1] 

Relative  velocity  at  tail 

^em  ( ^ ) 

[m  s-1] 

ECM  water  velocity 

4>em 

V 
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Figure  7.  Vertical  wavenumbei  spectra  of  measured  data  from  28  PATC1IEX  drops, 
4-8  MPa  (wavenumbers  assume  0.01  MPa  =  1  in).  Spectra  are  for  total  (east+north) 
variance;  indicated  noise  levels  are  also  for  combined  components.  Left:  ECM  velocity 
spectrum  (light  line)  has  spikes  near  A'rot  and  2A’ri>t,  and  closely  follows  the  GM76 
spectrum  (dashed  line)  computed  with  jV  =  0.00343  s-1;  ACM  spectrum  (heavy  line) 
has  a  deep  notch  near  0.2  cpm.  Right:  accelerometer  spectrum  peaks  near  0.11  cpm. 
The  variance  removed  by  digital  low-pass  filtering  of  the  data  has  been  recovered  in 
these  spectra  (see  Appendix  A). 


^Uido  ,_s  ZUJ  /  e<J> 


UNIVERSITY  OF  WASHINGTON  •  APPLIED  PHYSICS  LABORATORY. 


IV.  ACM  RESPONSE  MODEL 

To  produce  oceanic  velocity  profiles  from  the  ACM  data,  we  must  know  how 
the  instrument  was  moving.  Our  response  model  considers  the  motion  as  the  sum  of 
the  vehicle’s  oscillation  about  its  center  of  gravity  and  its  gross  movement  as  a  point 
mass  (at  CG). 

The  first  step  of  processing  removes  tilt  from  the  ACM  data,  producing  the  hor¬ 
izontal  relative  velocity  that  would  have  been  measured  had  MSP  remained  perfectly 
upright.  The  contaminating  signals  include  the  component  of  the  fall  rate  detected 
by  the  tilted  ACM  and  the  velocity  of  the  ACM  as  tilt  oscillations  sweep  it  through 
the  water.  Two  model  parameters  are  introduced  that  allow  the  necessary  corrections 
to  be  inferred  from  the  accelerometer  data. 

The  second  processing  step  uses  the  corrected  relative  velocities  to  estimate 
horizontal  forces  on  the  vehicle.  Integration  of  the  induced  accelerations  yields  MSP’s 
horizontal  velocity,  which  is  then  added  to  the  relative  velocity  to  form  the  oceanic 
velocity.  The  model  utilizes  one  force  near  the  nose  and  another  near  the  tail;  both 
vary  linearly  with  the  relative  velocity.  Each  force  is  specified  by  three  parameters: 
its  factor  of  proportionality,  its  location,  and  the  size  of  its  smoothing  window.  The 
forces  are  computed  at  their  respective  locations  but  are  then  applied  in  combination 
at  the  vehicle’s  center  of  gravity. 

A.  Vehicle  Kinematics 

To  consider  motion  confined  to  the  arz-plane,  let  x  and  z  be  the  horizontal  and 
vertical  unit  vectors  (Figure  8).  The  ocean  velocity  is  defined  as 

Vwi(i)  =  Ww«l(s)x  +  W'wat(*)£  (1) 


and  is  assumed  not  to  vary  during  an  MSP  drop.  The  vehicle’s  center  of  gravity,  P, 
has  location 


X,nsp(0  =  xp(t)x  +  zp(t)z 


(2) 


and  velocity 


( 0  —  ^  Xnisp(t)  Vnisp(f )  X  +  M  rosp(f)  Z 


(3) 


at  time  t.  The  fall  rate,  Wnisp  —  dzp(t)/df ,  is  positive  downward. 

When  the  vehicle’s  longitudinal  axis  (centerline)  is  rotated  anticlockwise  from 
vertical,  the  tilt  angle,  &(t),  is  positive.  An  additional  component  of  vehicle  motion 
occurs  as  the  tilt  changes.  In  particular,  the  velocity  of  the  ACM  is 


dO  . 

^msp(C  T  Lv  n  , 


(4) 
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Figure  8.  Definition  sketch  for  MSP  position  and  motion  at  time  1,  with  ACM  at 
depth  z(t).  Tilt  position,  0((),  and  oscillation,  dO/dt ,  are  shown  with  positive  sense. 
Unit  vectors  x,  z,  and  h(6)  point  in  the  positive  horizontal,  vertical,  and  transverse 
directions.  Point-mass  reference  P,  located  at  Xmsp (0  =  xp(l)x  +  zp(t)z,  moves  with 
velocity  v^pf/)  =  t'msP(i)x  +  WmspCO*-  At  distance  Lv  >  0  down  from  P,  the  ACM 
moves  with  an  additional  Lv  dO/dt  n  due  to  the  oscillating  tilt. 


where  Lv  is  the  distance  from  the  ACM  to  P  and 

fi(0)  =  cosO  x  —  sin#  z 


(5) 


is  the  transverse  unit  vector. 

The  ACM  detects  the  relative  motion,  along  the  direction  of  n,  of  the  water  at 
depth 

z{t)  =  zp{t)  +  Lvcosd  .  (6) 
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In  addition  to  the  horizontal  component  of  the  relative  velocity. 

t>rei .(0  =  vw«t(*(<))  -  *w(0  »  CO 

the  tilted  ACM  also  senses  part  of  the  vertical  component, 

HUO  =  W-tWO)  -  »W(f)  •  (8) 

Along  with  the  contribution  from  oscillations,  and  with  the  small-angle  approxima¬ 
tions  sin  0  =  0  and  cos  0  =  1,  the  total  ACM  measurement  is 

t’«cm (0  =  Vro,(<)  -  WU0*(0  -  (9) 

horiz  tilt  oscil 

neglecting  any  flow  distortion  at  the  transducers.  It  is  the  first  term,  the  relative 
horizontal  flow,  that  we  retain  for  the  response  model,  and  so  we  need  to  estimate 
the  tilt  and  oscillation  terms  to  remove  them  from  the  measured  ACM  data. 

We  infer  vehicle  tilt  from  the  accelerometer  data.  With  La  as  the  distance  from 
P  to  the  sensor  and  gravitation  acceleration  <7  =  9.8  m  s-2,  the  combined  signal  is 

,  ,  ,  dWnlsv  „  dvmsn  ....  dO  d26 

flmea»(0  ~  gO  ~  - ^ - &  +  -  ~  (Wnup  +  T’msp  9)  —  +  £a  09) 

tilt  horiz  frame  centrip 

assuming  small  angles.  The  first  term  is  the  vertical  acceleration  detected  by  a  tilted 
sensor  (note  that  <7  dW^p/dt).  In  Appendix  E  we  show  that  it  dominates  the 
remaining  terms  in  (10):  the  lateral  acceleration  of  MSP,  the  acceleration  perceived 
as  the  sensor  ;  reference  frame  changes  orientation  with  the  oscillating  tilt,  and  the 
centripetal  acceleration  of  the  sensor  relative  to  P. 


B.  Vehicle  Dynamics 

The  motion  of  the  water  past  MSP  leads  to  various  forces  that  push  the  vehicle 
horizontally  and  vertically  and  cause  it  to  tip.  These  forces  are  discussed  while 
developing  expressions  for  their  balances  in  Appendix  C.  Here,  to  set  up  the  response 
model,  we  summarize  the  horizontal  balance. 

As  illustrated  in  Figure  9,  the  nose  and  tail  forces, 

Fn  =  irr2NpVN  |Wrel|  and  Fr  =  vr(r2  -  r2N)pVT  |W'«,|  ,  [N]  (11) 
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Figure  9.  Definition  diagram  for  vehicle  dynamics  and  response  model  parameters. 
With  MSP  constrained  to  remain  upright,  the  entire  vehicle  moves  with  horizontal 
velocity  Origin  of  the  upward-pointing  .<s  axis  corresponds  to  reference  depth 

z  =  z(t)  of  the  ACM.  Variations  along  MSP  of  horizontal  water  velocity,  uw;,t(2-s), 
determine  those  in  relative  velocity,  uW4t(z-5)-umsp(f)  (shaded).  Defined  specifically  at 
the  ACM  (s=0)  is  the  relative  velocity  vrei(t).  Model  force  F/v  or  F-j  is  proportional  to 
relative  velocity  V/v(f)  or  Vj -(t),  the  average  over  a  window  of  size  or  Dr  centered 
Zn  or  Z,j  above  the  ACM.  Distances  Lv  and  La  separate  the  ACM  and  accelerometers 
from  the  reference  P. 
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pusli  MSP  horizontally.  The  total  vehicle  mass,  M,  includes  the  water  trapped  within 
the  outer  skin,  since  it  moves  along  with  the  solid  structure.  Furthermore,  MSP 
transmits  its  acceleration  to  some  volume,  Te,  of  the  exterior,  adjacent  water;  this 
retarding  action  is  known  as  the  added  mass  effect  (see  Carson  and  Simpson,  1978). 
Altogether,  the  balance  of  forces  yields  the  horizontal  acceleration  equation 

^  =  {f„  +  ft)/(m  +  p  re).  (12) 

Parameterization  and  integration  of  (12)  is  the  crux  of  the  response  model. 

C.  Model  Formulation 

The  primary  task  of  the  model  is  to  generate  vnisp(t),  the  instrument’s  horizontal 
velocity,  which  combined  with  the  ACM  measurements  produces  the  water  velocity. 
Because  the  observed  tipping  and  wobbling  arc  of  low  amplitude,  they  should  not 
appreciably  affect  the  gross  horizontal  motion  of  the  vehicle.  Therefore,  we  simplify 
the  dynamics  by  assuming  that  MSP  remains  upright  as  it  responds  to  the  water 
flowing  past.  However,  tilt  and  oscillation  effects  must  be  removed  from  the  ACM 
data,  j Vm(0'  supply  horizontal  relative  velocities  to  the  model.  These  corrections 
come  from  the  accelerometer  data,  amclUi(f).  The  two  horizontal  velocity  components 
(east  and  north)  result  from  separate  and  independent  runs  of  the  model. 

We  refer  to  the  position  of  MSP  at  time  t  by  the  depth  z  ■=  z(t)  of  the  ACM,  with 
z  positive  downward.  Variables  related  to  the  instrument’s  motion  or  measurements— 
those  defined  as  functions  off  in  Table  5  —are  henceforth  expressed  as  functions  of  z(t) 
or,  equivalently,  of  z.  For  instance,  uIlls,,(r)  =  vlus, ,(.?(<))  =  *Wp(0  is  the  horizontal 
velocity  of  MSP  when  the  ACM  is  at  depth  z  —  z{t).  Note  that  our  data  are  digitized 
on  1-kPa  pressure  grids,  nearly  equivalent  to  0.1- m  depth  grids. 

To  compute  vn.\(z)  from  the  ACM  data,  we  need  to  know  the  relative  vertical 
flow  and  the  vehicle  tilt.  Given  that  tilt  is  dominant  in  (10)  and  that  g  ^  dWmsp/ dt, 
we  set 

d[z)  =  aluaa{z)/g  [radians],  (13) 

We  estimate  the  fall  rate  from  the  pressure  ora  and  neglect  VFwai (z)  to  set  W/',.e| ( z )  = 
—  U'„isp(~ )  from  (8).  Now,  we  rearrange  (9;  <tud  compute 

Wrel(^)  =  Vicn.fz)  -  BjWuap{z)0{z)  +  Lv  [ill  S-1]  (14) 

from  the  0.1-m  ACM  data.  Parameter  if/  is  introduced  to  allow  investigation  of  the 
tilt  (fall-rate  bias)  correction.  Hayes  et  al.  (1984)  and  Schmitt  et  al.  (1988)  ignore  this 
term  and  thus  implicitly  have  Bj  —  0,  whereas  our  development  sets  Bj=  1;  we  alter 
Bj  between  these  extremes  in  tuning  and  in  assessing  its  effect  on  the  model.  Removal 
of  oscillation  from  the  ACM  measurements  via  (14)  is  vital  to  the  ultimate  success 
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of  the  vehicle  response  model.  Fortunately,  the  final  model  results  are  insensitive  to 
the  0. 1-0.2  m  uncertainty  in  Lv  arising  from  our  imprecise  knowledge  of  the  location 
of  P. 

Supplied  with  the  horizontal  relative  velocity  at  the  ACM,  we  need  only  t,„isp(c) 
to  determine  vwa,t(z)  from  (7).  As  shown  in  Figure  9,  we  attach  the  origin  of  the 
upward  s-axis  to  z  =  z(t),  so  that  the  relative  velocity  across  MSP  at  distance  s  above 
the  ACM  is  i Vat ( c  —  s)  —  vmsp(*),  which  for  a  =  0  is  urei(z).  We  define  the  time  and 
depth  of  the  drop’s  first  measurement  as  t  =  0  and  z(t)  =  z0  .  Our  initial  conditions 
are  vwaL(r  <  z0)  =  0  and  fmsp^o)  =  0.  Input  erei(z)  and  output  records  are  on  grids 
of  As  =  0.1  m. 


Our  neglect  of  VFw»t  implies  |M'rei(*)|  =  W7rasp(z)  =  dz/dt,  which  allows  us  to 
simplify  the  horizontal  acceleration  equation,  (12),  into 


dUmsp  (  Z  ) 


=  HNVN(z)  +  HrVr(z) 


with  force  parameters 


H\<  = 


nr 


N 


(Af/p+U) 


and  Hr  = 


7r(rf 


N 


(M/p  +  re) 


[nT1] 


and  smoothed  relative  velocities 


II 

2 

1 

Djv  ■ 

yZ/,'+0.5l)jv 

/  (  ^  Vat  (  ~  $ )  T’jnsp  ( Z )  )  ds 

J  Za,--0.5I)jv 

[m  s  1 

I't(z)  — 

1 

Dj  . 

yZ/'+O  .bDr 

/  (  tVat(-  -  S)  -  Voupiz)  )  ds 

IZt-O.SVt 

[m  s"1] 

(15) 

(16) 


(17) 


Parameters  Zk  and  Zj  specify  the  separation  between  the  ACM  and  the  centers  of 
application  of  Fs<  and  Fjy  respectively  (Figure  9),  and  Dpj  and  Dj  are  the  lengths  of 
the  corresponding  smoothing  windows. 

In  Appendix  D,  we  discuss  the  numerical  method  by  which  we  update  v,nsP(2) 
via  (15)  and  (17)  with  previously  computed  (or  initialized)  values  of  vvaX ('— .s)  and  the 
last  vehicle  velocity,  vtusp(z~  Ac) .  Given  the  relative  velocity  from  the  tilt-corrected 
ACM  and  the  vehicle  velocity  from  the  model,  we  can  finally  compute  the  water 
velocity 

^wal(^)  —  ^msp('<')  T  Urei(ir)  .  (18) 


The  initialization  uwat(c<r0)  =  0  establishes  a  depth-independent  offset  in  the 
oceanic  velocity  profile,  but  it  can  be  removed  with  knowledge  of  the  absolute  velocity 
at  some  depth  (from  a  nearby  mooring,  for  instance).  The  initial  lack  of  information 
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also  leads  to  a  startup  transient  in  (17),  which  vanishes  after  MSP  has  descended  a 
few  lengths.  Also,  this  portion  of  a  drop  generally  suffers  from  surface  wave  effects 
and  fall-rate  spinup. 

D.  Parameter  Values 

As  shown  in  Table  6,  we  have  assigned  nominal  values  to  the  eight  model  pa¬ 
rameters.  The  nose  and  tail  parameters  Zn,  D jv,  Zr ,  and  Dj  are  estimated  where 
the  vehicle’s  radius  changes  (as  in  Figure  9).  We  compute  H n  and  Hj  by  inserting 
known  or  estimated  dimensions  in  (16).  We  get  Lv  by  placing  P  near  MSP’s  center  of 
mass  and  set  Bf  =  1  consistent  with  our  formulation  in  (9).  Running  the  model  with 
the  nominal  values  leads  to  unsatisfactory  spectral  results,  however,  and  so  we  must 
determine  three  values  empirically:  Zr  to  match  the  0.2-cpm  notch,  Hh  to  attain 
the  low  wavenumber  shear  established  by  the  ECM,  and  Dr  to  dampen  the  vehicle 
response  at  high  wavenumbers.  We  can  attain  further,  more  subtle  improvements  to 
the  model  results  by  fine  tuning  the  remaining  parameters. 

To  find  optimal  parameter  values  for  a  given  MSP  configuration,  we  could  run  the 
model  using  various  sets  of  values  and  compare  the  spectra  of  the  resulting  profiles.  A 
quicker  and  more  powerful  approach  involves  analytic  transfer  functions  derived  from 
Fourier  transforms  of  model  equations  (in  Appendix  E).  Those  given  in  Table  5  are 
defined  for  the  response  of  each  quantity  to  the  water  velocity,  vWAl  Transfer  functions 
involving  only  the  modeled  point-mass  response,  A(Jb),  are  excellent  predictors  of  the 
corresponding  results  of  the  computational  model.  This  is  not  true  if  the  tilt  response, 
B(k),  is  involved,  since  it  is  derived  from  the  modeled  tilt  (Appendix  C),  which  differs 
significantly  (especially  for  k  >  0.4  cpm)  from  the  observed  tilt  as  inferred  from  (13). 

Recall  the  two  main  steps  of  our  model:  First,  generate  by  correcting 
for  tilt  and  oscillation;  second,  use  to  simulate  i’nup,  and  combine  them  to  produce 
uWat •  Because  the  first  step  involves  Lilt,  we  must  choose  values  for  Bj  and  Lv  by  trial 
and  error.  The  second  step  is  independent  of  tilt,  and  so  we  tune  the  six  potential 
force  parameters  using  CT{k),  the  transfer  function  for  the  corrected  relative  velocity 
Urci.  Specifically,  given  from  the  relative  velocity  profiles,  we  compute  CT{k)  for 
a  choice  of  parameter  values  and  generate 

=  *r7ia(*)i2 

for  a  trial  oceanic  spectrum.  The  result  closely  represents  the  spectrum  of  owat  profiles 
that  would  be  generated  from  vn\  via  (15),  (17),  and  (18). 

Figure  10  compares  v^,  for  drop  0096  with  as  computed  from  (14)  with 
nominal  values  of  Bj  and  Lv  (Table  6).  The  tilt  and  oscillation  terms  have  amplitudes 
of  1  and  2  mm  s-1.  Although  small,  these  corrections  are  significant  at  scales  near 
the  instrument  length.  Figure  11  shows  their  effect  on  the  measured  PATCI1EX 
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Table  6.  Response  model  parameters.  Nominal  values  are  those  expected  from  mea¬ 
sured  or  estimated  dimensions.  Optimal  values,  with  which  the  model  is  run,  are  de¬ 
termined  with  help  from  the  transfer  functions.  Values  in  parentheses  are  for  the  small 
tail. 


Parameter 

Description 

Nominal 

Value 

Optimal  Value 

Large  Tail  Small  Tail 

Zk 

M 

Distance  to  nose  force 

0.50 

0.50 

0.50 

ZT 

[m] 

Distance  to  tail  force 

3.40 

3.98 

4.00 

D  K 

M 

Nose  window  size 

0.60 

0.80 

0.80 

Dr 

[m] 

Tail  window  size 

0.80 

1.40 

1.60 

II N 

[nr1] 

Nose  force  factor 

0.18 

0.28 

0.32 

Hr 

[nr1] 

Tail  force  factor 

1.14  (0.54) 

1.16 

0.52 

Lv 

[nr] 

Distance  from  ACM  to  P 

1.9 

1.8 

2.0 

Dj 

Fall -rate  bias  factor 

1.0 

0.9 

0.6 

A 

[s'1] 

Tilt  damping 

1.18  (0.84) 

0.6 

1.0 

00 

[r2] 

Tilt  restoring 

0.70  (0.75) 

1.1 

1.4 

Qn 

[m-1  a”1] 

Nose  tilt  factor 

0.07  (0.09) 

0.05 

0.065 

Qt 

[in-1  s-1] 

Tail  till  factor 

0.55  (0.35) 

0.95 

0.95 

Ltt 

[m] 

Distance  from  acce).  to  P 

1.2 

spectrum,  4>'a‘a".  lu  the  corrected  spectrum,  4>,cl,  the  0.2-cpm  notch  is  broader  and 
half  as  deep,  amplitudes  of  the  other  peaks  and  valleys  are  diminished,  and  spectral 
levels  at  k  >  0.7  cpm  are  lowered.  Oscillation  accounts  for  most  of  the  difference, 
although  tilt  contributes  at  k  <  0.2  cpm.  These  basic  results  occur  for  all  MSP 
datasets. 

Id  complete  our  transformation  from  measured  to  oceanic  velocities,  we  must 
find  values  for  the  six  parameters  of  the  point-mass  model.  In  Appendix  E,  we  derive 
the  expression 


Hre~ikZr 

Cr{k)  =  1  -  77-1-77— nr  sinc(*D,v/2)  -  -  -j-ry—rr  sinc(fc£>r/2)  (19) 

i/yv  "4“  Ti  7  -j-  IK  i/yv  Il'j'  -{-IK 


where  k  is  the  wavenumber  in  radians  per  meter  and  sincx  =  x_1  sin  x  .  Starting  with 
t  he  nominal  values,  we  generate,  oceanic  spectra  from  C\  and  and  iteratively  tune 
the  parameters  until  the  result  satisfies  our  slated  requirements:  a  smooth  spectrum, 
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Figure  10.  Components  of  ACM-measured  velocity,  PATCIIEX  drop  0096:  (light 

line)  is  composed  of  (dark  line),  tilt  ( BjWnay>0 ),  and  oscillation  (-L^dO/dl).  For 
display,  tilt  and  oscillation  terms  were  multiplied  by  5  and  ofTset  by  0.04  and  0.06  m  s-1. 
Aside  from  a  few  instances  (e.g.,  near  6.75  MPa  for  the  east  components),  v^x  and  nrcj 
are  indistinguishable. 


devoid  of  obvious  instrumental  effects,  which  is  close  to  ECM  levels  at  low  wavenum¬ 
bers.  In  Figure  12  we  present  the  optimally  tuned  versions  of  Cr  for  our  two  tail 
configurations,  and  in  Table  6  we  list  the  corresponding  sets  of  parameter  values.  We 
truncate  the  display  at  k  —  1  cpm,  since  CT(k)  1  at  higher  wavenumbers.  Unless 
otherwise  noted,  the  remaining  discussion  pertains  to  the  configuration  used  through 
Tropic  Heat  2  (large  tail). 

Figure  13  compares  the  Cr  computed  from  (19)  when  using  only  the  tail  force 
(f/yv  =  0)  with  that  computed  when  using  only  the  nose  force  (Hr  =  0).  The 
dominance  of  the  tail  force  is  obvious,  judging  by  the  similarity  of  CT  to  that  of 
the  full  model.  Its  position  Zj  fixes  the  wavenumbers  of  the  notch  and  subsequent 
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Figure  11.  Shear  spectra  (at  intermediate  stage)  of  28  PATCUEX  drops,  4-8  MPa, 
computed  as  —  (2trfc)2$v.  Spectra  for  measured  ACM  records  are  shown  by  the 
light  line  and  tilt-corrected  spectra  by  the  heavy  line  (shaded  at  95%  confidence  limits). 
Also  shown  are  spectra  of  the  tilt  and  oscillation  corrections  and  of  the  GM76  model. 
Spectra  are  truncated  at  1  cpm,  beyond  which  vehicle  motion  has  little  effect. 
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Figure  12.  Transfer  functions  Cr(k)  from  ocean  velocity  to  relative  horizontal  velocity 
at  ACM.  as  computed  with  the  optimal  parameter  values  in  Table  6:  Upper  panel  is 
squared  amplitude;  lower  is  phase,  between  0.1  and  0.3  cpm,  the  response  for  the  large- 
tail  array  used  up  through  Tropic  Heat  2  (solid  line)  is  more  severe  than  that  for  the 
small-tail  array  used  after  Tropic  Heat  2  (dashed  line). 


spectral  features,  and  its  strength  Hr  sets  their  amplitudes.  Its  window  size  Dr 
attenuates  high-wavenumber  features  by  reducing  the  response  of  the  vehicle  to  water 
motions  of  small  scale.  In  regard  to  the  nose  force,  its  strength  Hu  establishes  the 
low- wavenumber  response,  its  position  Z ^  plays  a  minor  role,  and  its  window  size 
Dk  has  almost,  no  effect. 

Evaluation  of  sensitivities  of  the  model  to  each  of  the  eight  parameters  is  a  by¬ 
product  of  tuning.  As  long  as  realistic  values  are  assigned  to  B;  and  Lv ,  the  success 
of  the  model  is  insensitive  to  their  exact  value,  Without  the  oscillation  correction 
implied  by  even  the  most  carefully  chosen  values  of  Zr  and  Hr  yield  that 
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wavenumber  /  cpm 

Figure  13.  Transfer  functions  Cr( k)  for  the  large-tail  array  computed  exclusively  for 
the  tail  force  (solid  lines)  by  setting  7//v  =  0  and  exclusively  for  the  nose  force  (dashed 
lines)  by  setting  II r  —  0.  With  only  the  tail  force,  CTr ( A* )  resembles  the  complete  two- 
force  version,  aside  from  the  difference  in  amplitude  for  k  <  0.1  cpnt.  With  only  the  nose 
force,  Cr(k)  loses  the  shape  of  the  complete  version,  although  the  amplitudes  converge 
at  k  <  0.04  cpm. 


are  choppy  at  k  >  0.2  cpm.  Given  the  corrected  <1*™',  we  can  focus  upon  the  six 
parameters  in  (19).  The  point-mass  model  is  most  sensitive  to  Zj,  coupled  with 
//■;•,  since  these  parameters  must  account  for  the  remaining  0.2-cpm  notch.  Because 
the  notch  is  deep  and  narrow,  its  position  must  be  closely  matched  by  Cr  for  the 
resulting  to  be  acceptably  smooth.  The  nominal  Zj  =  3.5  m  fails  the  test,  as 
do  values  more  than  0.05  m  away  from  the  optimal  value  of  3.98  m  (which  places  Fy 
up  in  the  wake  of  the  brushes).  Sensitivity  to  the  other  parameters  is  more  subtle, 
since  variations  in  their  values  affect  broad  spectral  levels  rather  than  shapes.  Two  of 
them,  Dj  and  //*>,  differ  noticeably  from  the  nominal  predictions.  The  enlarged  tail 
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window  is  needed  to  squash  the  response  beyond  0.4  cpm;  the  stronger  nose  force  is 
needed  to  bring  up  toward  the  level  of  the  ECM  at  low  wavenumbers.  Finally, 
the  position  of  F /v  slightly  affects  the  notch  depth. 

Minor  variations  in  the  optimal  values  of  the  parameters,  other  than  Zj,  do  not 
significantly  alter  Cr  (or  Cm).  In  Appendix  E,  we  fine  tune  Lv  and  Bj  using  the 
transfer  function  Hm(k)  from  14^  to  anie*».  We  compensate  for  any  resulting  effect 
on  4>'sel  by  tweaking  Hr  or  Zyv . 

Wc  summarize  the  action  of  MSP  as  it  filters  the  water  motion,  comparing  its 
point-mass  response  and  its  tilt  response.  The  large-scale  motions  carry  the  vehicle 
as  a  point  mass,  with  the  nose  force  dominating  for  k  <  0.04  cpm.  From  0.05  to 
0.15  cpm,  the  primary  effect  comes  from  the  tail  force  and  its  separation  from  the 
ACM.  Tilt  and  oscillation  start  to  contribute  at  k  >  0.1  cpm;  between  0.2  and 
0.4  cpm,  the  effect  of  oscillation  reinforces,  and  nearly  equals,  that  of  the  tail  force. 
Since  the  point-mass  response  is  weak  for  k  >  0.4  cpm,  the  oscillating  tilt  of  the  ACM 
dictates  the  difference  between  v ^  and  uwat  at  the  smallest  scales. 

After  Tropic  Heat  2,  we  designed  the  smaller  tail  to  lower  II j  via  (16).  As 
predicted  by  (19)  and  reflected  in  Figure  12.  the  weaker  tail  force  led  to  a  broader 
and  shallower  notch  in  Relative  to  that  for  the  large  tail,  Cr  for  the  small-tail 
configuration  is  less  sensitive  to  variations  in  the  tail  parameters. 
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V.  RESULTS 

A.  Velocity  Profiles  and  Shear  Spectra 

Ocean  velocity  profiles  from  the  ACM  response  model  nicely  track  the  ECM 
results,  aside  from  some  large-scale  divergence.  Figure  14  demonstrates  the  model's 
success  for  PATCHEX  drop  0096  and  also  for  a  Florida  Straits  drop  containing  consid¬ 
erable  mean  shear  and  strong  fine-scale  shear.  As  anticipated,  the  agreement  between 
profiles  degrades  at  scales  smaller  than  the  rotation  length. 

Spectral  results  are  a  critical  test  of  the  model's  performance.  After  running 
the  28  PATCHEX  drops  through  the  model,  we  form  the  shear  spectrum  of  the  nWi4t 
profiles  for  comparison  with  the  ECM  spectrum  (Figure  15).  At  low  wavenumbers, 
up  to  0.04  cpm,  the  spectra  are  similar  in  level  (both  being  close  to  GM76)  and  are 
highly  coherent  with  near-zero  phase.  Around  AVot  and  2  A'rol,  coherence  is  minimal 
and  spikes  in  4>'“‘  reach  twice  the  ACM  levels.  At  higher  wavenumbers,  the  spectra 
remain  fairly  coherent  and  close  in  phase  despite  the  higher  ECM  levels  attributed  to 
rotation-related  contamination.  Coherence  is  lost  at  k  >  0.5  cpm  as  noise  dominates 
the  ECM.  The  ACM  spectrum  is  reasonably  smooth  between  0.1  and  1  cpm,  aside 
from  minor  irregularities  at  0.4-0. 5  cpm.  To  emphasize  the  importance  of  the  tilt 
effects,  we  include  the  spectrum  of  profiles  created  without  correcting  for  tilt  and 
oscillation,  that  is,  by  running  the  model  with  Bj  =  Lv  —  0  and  with  Zj  and  Hj 
chosen  so  that  the  tail  force  alone  compensates  for  the  diep  notch  in  $Jcnl.  Now,  the 
spectrum  is  excessively  choppy  and  uneven  in  the  high-wavenumber  band,  with  the 
coherence  and  phase  relation  to  the  ECM  breaking  down  around  0.2  cpm. 

Many  features  of  the  PATCHEX  results  also  appear  in  the  more  energetic  Florida 
Straits  data,  taken  with  the  small-tail  configuration.  The  ACM  and  ECM  spectra  are 
coherent  with  comparable  levels  at  low  wavenumbers;  their  coherence  fades  near 
and  2A'rot,  rises  at  higher  wavenumbers,  and  vanishes  as  noise  overtakes  the  ECM 
past  0.4  cpm  (Figure  16).  Unlike  for  the  large  tail,  oscillation  correction  has  little 
effect  on  coherence  near  0.2  cpm,  although  it  remains  a  vital  step  toward  generating 
a  smooth  ACM  spectrum. 


We  check  for  consistency  in  the  ACM  and  airfoil  levels  in  the  band  between 

1er>  *-»-»  lull  i  pf ri  i  1  ('  rv  /-k/»  t  r U  *-»■»  «  »-rtl  >  -i  L  1  /->  •%  rl  O  4  /-,  Off  A. <irU  wn  A  f"*  A  <4 

v.iiwv  tiAiivsii  uuvvmt  tuiauic,  a.iiU  tu  L\)  CjJUi,  wiiciC  ia«  /VvaVI 

data  reach  their  noise  floor.  For  the  most  direct  comparison,  we  utilize  data  from 
the  ACM  axis  that  is  more  closely  aligned  with  the  airfoil  probes.  The  spectra 
generally  agree  to  within  a  factor  of  1.5,  although  there  is  a  tendency  for  airfoil 
levels  to  exceed  ACM  levels  by  a  ratio  that  increases  with  greater  microscale  shear 
variance  or  higher  temperature  (Appendix  B).  When  we  superimpose  airfoil  spectra 
onto  processed  ACM  spectra,  O**1,  we  anticipate  some  deviation  since,  in  the  latter, 
subsampling  leads  to  aliasing  at  k  >  3.5  cpm  and  correcting  for  oscillation  reduces 
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Figure  14.  Water  velocity  profiles  from  ACM  response  model  (heavy  line)  and  ECM 
results  (light  line)  for  (a)  PATCHEX  drop  0096  and  (b)  Florida  Straits  drop  0241.  All 
profiles  are  ofTset  for  comparison,  with  those  in  (b)  forced  to  be  near  zero  at  7  MPa. 
For  the  east  component  at  superrotation  scales,  the  ACM-ECM  rms  difference  in  (a) 
is  7  mm  s-1;  in  (b),  it  is  15  mm  s_1  from  1  to  6.4  MPa  and  7  mm  s-1  between  3 
and  6  MPa.  For  the  north  component,  the  difference  in  (a)  is  4  mm  s-1;  in  (b),  it  is 
19  mm  s-1  from  1  to  6.4  MPa  and  13  mm  s-1  between  3  and  6  MPa.  Note  that  the 
velocity  scale  in  (a)  is  a  tenth  of  that  in  (b). 


levels  at  k  <  2  cpm.  Nonetheless,  there  is  usually  a  crossing  or  close  matching 
somewhere  between  1  and  3.5  cpm  (Figure  17). 

B.  Discussion 

The  results  demonstrate  our  success  in  combining  data  from  the  ECM,  ACM,  and 
airfoil  probes  to  resolve  velocity  and  shear  in  the  band  from  0.01  to  100  cpm.  The 
response  model  generates  shear  spectra  that  are  consistent  with  the  ECM  spectra 
at  low  wavenumbers  and  with  the  airfoil  probe  spectra  at  high  wavenumbers  and 
reasonably  accounts  for  features  of  measured  ACM  spectra  in  the  intermediate  band. 
However,  there  are  uncertainties  that  require  comment. 
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Figure  14,  cont. 


While  the  agreement  of  ACM  and  ECM  results  is  satisfying,  we  cannot  consider 
the  ECM  an  absolute  standard,  and  so  seek  further  validation.  Peters  et  al.  (1991) 
compare  ACM  profiles  from  several  MSP  drops  through  the  Equatorial  Undercurrent 
(EUC)  with  simultaneous  measurements  by  an  RDI  acoustic  Doppler  current  profiler 
(ADCP)  (the  ECM  was  ineffective  owing  to  slow  MSP  rotation  and  the  proximity  of 
the  magnetic  equator).  Considering  mean  profiles  beneath  the  EUC  core,  the  shear  is 
10%  higher  in  the  ACM  data  than  in  the  ADCP  data  (which  was  prone  to  systematic 
errors  in  this  regime).  In  contrast,  ACM  levels  for  PATCHEX  are  about  10%  lower 
than  the  ECM  levels  (Figure  15).  We  feel  that  our  low- wavenumber  tuning  (via  Hn) 
is  a  good  compromise.  With  the  additional  concerns  from  our  use  of  nominal  ACM 
calibrations  and  neglect  of  lateral  drag  and  vertical  water  motions,  we  estimate  a 
combined  uncertainty  of  10%-15%  in  the  large-scale  ACM  results. 

At  scales  around  the  vehicle’s  length,  there  is  no  standard  for  objectively  judging 
the  ACM  profiles.  Between  0. 1  and  1  cpm,  we  can  only  require  that  spectra  are  smooth 
in  shape,  that  is,  that  the  model  has  decently  corrected  the  effects  reflected  in  the 
measured  spectra.  The  model  meets  this  criterion  fairly  well.  The  similarity  to  GM76 
spectra  and  coherence  with  the  ECM  spectra  are  also  encouraging.  Irregularities  that 
sometimes  appear  near  0.2  and  0.4-0. 5  cpm  might  imply  a  narrowband  uncertainty 
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wavenumber  /  cpm 

Figure  15.  Total  shear  spectra  of  28  PATCIfEX  diops,  4-8  MPa,  computed  from 
ocean  velocity  profiles.  Top  panel:  ECM  spectrum  (light  line,  truncated  at  0.4  cpm), 
standard  ACM  spectrum  4>*al  (heavy  line,  with  95%  confidence  limits),  tiltless  ACM 
spectrum  (light  line)  calculated  using  Bj  =  Lv  =  0  (see  text),  and  GM76  spectrum 
(dashed  line).  Bottom  panels:  Squared-coherence  and  phase  between  ACM  and  ECM 
velocities  (east  components),  95%  confidence  intervals  shaded.  Heavy  lines  are  for  stan¬ 
dard  ACM;  light  lines  (near  0.2  cpm)  are  for  ACM  processing  without  tilt  and  oscillation 
corrections.  The  phase  truncates  at  0.6  cpm,  where  the  coherence  vanishes. 
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Figure  16.  Total  shear  spectra  of  14  midchannel  Florida  Straits  drops,  2.4-6.45  MPa, 
computed  from  ocean  velocity  profiles.  Spectral  processing  is  similar  to  that  for  PATCH- 
EX,  except  that  transforms  are  taken  on  segments  of  2.7  MPa  rather  than  1.0  MPa.  Top 
panel:  ECM  spectrum  $*m  (light  line),  ACM  spectrum  (heavy  line),  and  GM76 
spectrum  (dashed  line,  A  =0.00504  s-1).  Bottom  panels:  Squared-coherence  and  phase 
between  ACM  and  ECM  (east  components). 
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Figure  17,  PATCH EX  North  \CM  ($7**)  and  airfoil  ($yuu'°)  shear-  spectra  computed 
from  5.75  to  9.25  Ml’a.  The  airfoil  spectrum  is  doubled  for  comparison  with  the  ACM 
total  shear  spectrum.  The  GM7G  internal  wave  model  shear  spectrum  (dashed  line, 
N  =  0.00315  s-1)  and  Nasmyth  turbulent  model  shear  spectrum  (dotted  line,  e  =  1.94  x 
10-!1  W  kg"1,  doubled  for  comparison)  are  plotted  for  reference.  The  ACM  spectrum, 
corrected  for  variance  removed  by  the  2-kPa  Bartlett  window,  tails  up  after  3.5  cpm 
owing  to  aliasing  (Appendix  A).  The  airfoil  spectrum  is  unreliable  before  1  cpm  and  is 
at  noise  levels  past  30  cpm. 


of  10% -30%,  based  on  their  deviation  from  a  smooth  spectrum.  However,  when  such 
spectral  features  are  figured  into  the  integrated  shear,  their  significance  is  diluted. 
Individual  profiles  can  suffer  from  effects  that  are  not  obvious  in  spectral  results. 
For  example,  the  model  occasionally  generates  5- in  oscillations  by  overamplifying 
features  from  the  measurements.  This  is  most  apparent  when  shear  at  neighboring 
scales  (2-10  in)  is  weak  over  some  depth  interval. 

For  high  wavenumbers,  we  consider  the  ACM  as  the  standard  since  it  is  more 
stable  than  the  airfoil  probes.  Their  spectral  levels  compare  favorably,  especially  when 
considering  raw  ACM  data  (Appendix  B).  The  ACM  measurements  are  affected  by 
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oscillation  but  not  by  point-mass  motion,  so  the  uncertainties  come  from  Lv  and 
from  calibration  and  noise  in  the  ACM  and  accelerometers.  We  are  satisfied  with  the 
microscale  results,  since  our  goal  is  to  patch  together  ACM  and  airfoil  spectra  toward 
resolving  shear  variance  out  to  the  viscous  rolloff. 

The  concerns  we  have  raised  thus  far  are  mainly  model  related.  Now  we  address 
uncertainties  introduced  as  the  ACM  measurements  are  transformed  into  geomag¬ 
netic  coordinates.  To  gauge  the  potential  effect  on  our  results,  we  perturb  zero 
offsets,  calibrations,  or  alignment  adjustments  and  reprocess  the  data.  For  tests 
run  on  PATCHEX  drops,  minor  perturbations  have  little  effect  on  4>*al.  Excessive 
modification  of  zero  offsets  (by  10  mm  s-1  in  contrast  to  expected  uncertainties  be¬ 
low  2  mm  s"1)  yields  the  strongest  response,  elevating  the  spectrum  tenfold  at  A’rot. 
Moderately  altering  the  sensitivity  (10%)  or  the  alignment  (10°)  of  a  single  ACM 
axis  increases  variance  by  5%-15%  in  the  band  from  0.17  to  0.25  cpm.  These  effects 
appear  predominantly  in  the  clockwise  component  of  rotary  spectra  (computed  in  the 
downward-looking  oceanographic  convention;  see  Gonella,  1972). 

C.  Noise  and  Errors 

The  ACM  noise  level  sets  a  lower  bound  on  the  noise  in  our  ocean  velocity 
profiles.  Unless  N  <  10~3  s"1,  spectra  conforming  to  GM76  stay  above  the  ACM 
noise  floor  (10~8  m2  s~2  cpm-1)  for  k  <  1  cpm.  Low-pass  filtering  the  raw  data 
reduces  the  rms  noise;  transforming  to  geomagnetic  coordinates  and  correcting  for 
oscillation  could  add  some  noise.  The  degree  of  uncertainty  in  the  model  results 
depends  on  the  scale.  For  scales  over  10-20  m,  where  velocity  variance  resides,  the 
effect  is  equivalent  to  a  10%-15%  calibration  uncertainty.  At  scales  exceeding  the 
rotation  length,  rms  differences  between  ECM  and  ACM  profiles  are  3-20  mm  s-1. 
For  2-  10  m,  the  vehicle  response  fluctuates  such  that  the  model  amplifies  some  scales 
and  attenuates  others.  Sensitivity  is  greatest  around  the  5-m  scale.  We  estimate 
uncertainties  of  10%  in  general  and  up  to  30%  at  spectral  irregularities.  To  assess 
uncertainty  in  a  specific  profile,  we  could  examine  the  effects  as  the  model  is  run  over 
a  range  of  parameter  values. 

Our  studies  rely  more  on  shear  and  Richardson  numbers  than  on  velocity  profiles. 
Noise  and  uncertainty  in  shear  depend  on  the  scale  over  which  it  is  computed.  If  larger 
scales  are  the  focus,  we  reduce  the  noise  by  filtering  the  velocity  data.  Regarding 
shear  spectra  noise  levels,  the  ACM  and  the  airfoil  probes  are  both  around  4  x 
10~7  s'2  cpm-1  at  1  cpm,  with  the  ACM  floor  sloping  up  thereafter  and  that  of  the 
airfoils  falling  slightly  out  to  30  cpm.  Uncertainty  in  shear  is  related  to  that  in  velocity; 
it  is  easier  to  quantify  at  large  scales  than  at  intermediate  scales.  Unfortunately,  shear 
is  concentrated  at  scales  where  the  response  model  is  most  questionable.  Again,  we 
could  make  specific  error  estimates  by  running  repeated  simulations  on  the  drops  in 
question. 
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Gregg  (1979)  presents  several  formulae  for  estimating  relative  uncertainty  and 
noise  in  computed  quantities.  The  relative  noise  a^j Ri  increases  with  noise  of  the 
temperature,  conductivity,  pressure,  and  velocity  sensors,  decreases  with  larger  inter¬ 
val  size,  and  increases  with  the  actual  value  of  Ri.  For  MSP,  noise  of  the  CTD  sensors 
contributes  much  more  than  that  of  the  ACM.  While  we  can  apply  the  relation  to 
estimate  <rRJ Ri  for  a  given  MSP  profile,  we  must  rely  on  repeated  simulations  to 
estimate  relative  uncertainty. 
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VI.  CONCLUSIONS  AND  SUMMARY 

High-quality  measurements  mark  the  sound  performance  of  MSP.  Success  comes 
from  the  combined  abilities  of  its  sensor  systems  coupled  with  its  advantageous  behav¬ 
ior  as  a  platform.  The  drag  brushes  slow  the  descent,  enhancing  resolution  of  temper¬ 
ature  inicrostructure.  With  a  fall  that  is  quiet  and  low  in  vibration,  MSP  can  detect 
a  wide  range  of  viscous  dissipation  rates.  Rapid  sampling  of  temperature  and  ve¬ 
locity  data  allows  broadband  comparisons  between  fine-scale  and  microscale  sensors. 
In  particular,  ACM  data  are  useful  in  assessing  the  veracity  of  airfoil  spectral  levels. 
The  ECM  provides  a  check  on  large  scale  ACM  results  without  deployment  of  extra 
equipment  (as  would  be  required,  for  instance,  with  acoustic  tracking).  Adjustable 
turning  blades  spin  the  vehicle,  meeting  a  requirement  of  the  ECM  method.  Rotation 
also  enables  determination  of  bias  and  misalignment  in  the  ACM,  accelerometer,  and 
magnetometer  channels.  A  single  set  of  accelerometers  is  sufficient  for  monitoring 
vehicle  till,  since  free  oscillation  is  suppressed  by  damping  and  off-resonant,  forcing 
at  the  fail.  However,  by  strengthening  lateral  forcing,  the  tail  array  of  brushes  and 
blades  adversely  affects  the  ACM  measurements. 

We  have  focused  this  report  on  our  treatment  of  the  ACM  data  because  the 
resulting  oceanic  profiles  and  spectra  are  central  to  MSP  investigations  of  shear  and 
mixing.  With  its  stabiliiv  and  low  noise,  the  ACM  can  resolve  motions  on  vertical 
scales  from  1  to  10  m  and  beyond.  It  measures  the  relative  How  past  its  transducers, 
and  thus  proper  interpretation  of  the  data  requires  knowledge  of  the  orientation  and 
motion  of  the  platform.  Determining  the  gross  vehicle  motion  is  the  most  difficult 
task,  for  which  we  rely  on  a  computational  model  based  on  our  understanding  of  the 
horizontal  dynamics  affecting  MSP.  At  large  scales,  forcing  at  the  nose  couples  MSP 
to  the  flow,  much  like  the  profilers  TOPS  and  HRP.  The  troublesome  aspect  of  the 
response  is  the  tendency  of  MSP,  because  of  enhanced  tail  forcing,  to  follow  motions 
of  scales  near  its  own  length.  This  behavior,  combined  with  the  swaying  of  the  ACM 
from  tilt  oscillations,  complicates  the  ACM  measurements  such  that  their  spectra  are 
deeply  notched  near  0.2  cpin  with  lesser  contortions  thereafter.  To  be  successful,  our 
model  must  account  for  such  features  in  the  measured  spectra..  While  ECM  results 
guide  us  at  scales  exceeding  the  rotation  length,  no  standard  exists  at  smaller  scales. 

'The  transformation  from  ACM  data  into  oceanic  velocity  profiles  occurs  in 
stages.  First,  we  use  magnetometer  data  to  resolve  the  ACM  and  accelerometer 
data  into  geomagnetic  coordinates.  Next,  we  remove  vehicle  till  and  oscillation  from 
the  ACM  measurements,  inferring  the  corrections  from  till -dominated  accelerometer 
data.  The  model  simulates  the  gross  horizontal  motion  of  the  vehicle  as  it  is  responds 
to  forces  at  the  nose  and  tail.  Three  parameters  specify  the  location  and  strength 
of  each  force  in  proportion  to  t he  oncoming  relative  flow.  We  manufacture  the  ve¬ 
hicle  velocity  by  integrating  the  modeled  accelerations  and  then  compute  the  ocean 
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velocity  by  adding  on  the  corrected  relative  measurements.  Rather  than  adhering 
to  supposed  or  theoretical  values  for  the  parameters,  we  adjust  them  empirically  by 
matching  analytic  transfer  functions  to  the  observed  response. 

Velocity  profiles  and  shear  spectra  produced  by  the  ACM  processing  satisfy 
most  of  our  stated  objectives  and  criteria.  At  large  scales,  the  ACM  results  agree 
well  with  ECM  profiles  and  spectra.  The  accelerometers  and  response  model  capture 
the  effects— oscillation  and  tail  forcing — that  shape  the  measured  spectra  between 
0.1  and  1  cpm.  Consequently,  our  ocean  shear  spectra  are  fairly  smooth  within  this 
difficult  band,  aside  from  minor  deviations  near  0.2  and  0.4-0. 5  cpm.  Airfoil  and  ACM 
spectra  usually  cross  or  overlap  between  1  and  3  cpm.  The  consistency  between  ECM, 
ACM,  and  airfoil  results  confirms  the  ability  of  MSP  to  resolve  oceanic  shear  spectra 
from  0.01  to  100  cpm. 

Overall,  we  are  satisfied  with  our  ACM  results,  but  we  do  have  some  comments 
and  concerns.  Considering  sensitivities  in  the  measurements  and  model  parameters, 
uncertainty  in  spectral  levels  is  10%  in  general,  but  is  potentially  higher  near  0.2  cpm 
and  near  the  rotation  wavenumber,  A'lot,.  Some  parameters  differ  significantly  from 
their  predicted  values:  a  strong  nose  force  and  an  upwardly  displaced  tail  force  are 
required  for  success.  In  the  spatial  domain,  some  portions  of  profiles  should  be  viewed 
with  suspicion.  Near  the  surface,  results  are  confounded  by  surface  waves,  vehicle 
spinup,  and  startup  transients  in  the  model.  Neglect  of  lateral  drag  on  the  main 
tube  may  cause  the  model  to  underestimate  acceleration  through  sharp  or  strong 
shear  features.  The  model’s  extreme  response  at  0.2  cpm  can  generate  unrealistic 
5-m  oscillations  in  velocity  or  shear  profiles  (more  of  a  concern  for  the  large-tail 
configuration  than  for  the  small-tail  one).  Finally,  velocity  profiles  contain  depth 
independent  offsets  and  extend  to  no  more  than  1200-m  depth. 

We  compare  MSI*  with  TOPS  and  HltP  in  terms  of  vehicle  behavior  and  treat¬ 
ment  of  ACM  data.  In  some  ways,  MSP  is  superior  in  acquiring  data:  a  slow,  quiet 
fall  improves  microstructurc  measurements,  rapid  sampling  allows  more  flexibility  in 
data  processing,  and  the  ECM  provides  au  onboard  check  on  ACM  results.  How¬ 
ever,  the  drag  brushes  and  turning  blades  that  control  the  descent  also  complicate 
the  ACM  measurements  by  causing  MSP  to  be  more  responsive  than  the  other  pro¬ 
filers.  Regarding  vehicle  tilt,  both  MSP  and  HRP  improve  upon  the  weakly  damped 
resonant  response  of  TOPS.  With  its  greater  length,  HRP  attains  a  strong  righting 
moment  and  large  moment  of  inertia,  leading  to  rapid  damping  of  its  free  oscillation. 
In  contrast,  the  tilt  response  of  MSP  is  more  forced  than  free,  as  opposing  nose  and 
tail  forces  tip  the  vehicle.  Also,  the  tail  effectively  stiffens  and  damps  the  resonant 
mode.  A  tailless  MSP  would  oscillate  badly,  judging  from  the  proximity  of  its  centers 
of  buoyancy  and  gravity.  The  TOPS  model  couples  the  equation  for  gross  vehicle 
motion  to  that  for  vehicle  tilt.  The  primary  use  of  its  noise-limited  accelerometer 
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measurements  is  in  tuning  the  tail  force  parameter.  Neither  MSP  nor  HRP  process¬ 
ing  includes  tilt  dynamics  in  modeling  vehicle  velocity.  Both  remove  oscillation  effects 
from  their  ACM  measurements,  but  only  the  MSP  corrects  for  the  component  of  fall 
rate  detected  with  tilt.  One  set  of  accelerometers  suffices  for  monitoring  MSP’s  tilt; 
HRP  has  sensors  at  both  ends  to  decouple  tilt  from  other  acceleration  signals.  All 
three  models  employ  a  linear  nose  force  and  ignore  lateral  drag,  vertical  water  motion, 
and  rotation  effects.  Peaks  near  the  pendulum  frequency  dominate  spectra  of  TOPS 
accelerometer  and  ACM  measurements;  it  is  only  in  the  transfer  function  between  the 
two  that  effects  of  tail  forcing  are  apparent.  The  HRP  model  omits  the  tail  force,  as 
oscillation  correction  seems  competent  to  account  for  spectral  irregularities.  The  MSP 
model  requires  a  tail  force  to  explain  the  0.2-cpm  notch  and  surrounding  features  in 
the  measured  ACM  spectra.  By  eliminating  dependence  on  fall  rate  and  combining 
several  physical  parameters,  we  employ  only  six  simple  parameters  in  our  vehicle  re¬ 
sponse  model.  We  use  model  transfer  functions  to  select  parameter  values,  to  assess 
model  sensitivities,  and  to  plan  modifications  toward  improving  vehicle  response. 

Much  of  the  design  of  MSP  involves  controlling  and  monitoring  vehicle  motion. 
The  array  of  drag  brushes  and  turning  blades  affects  all  aspects  of  the  response — 
vertical,  horizontal,  and  tilting.  The  large-tail  configuration  met  our  goal  of  a  slow 
descent  but  adversely  affected  the  ACM  measurements  by  making  MSP  respond  too 
quickly  to  horizontal  flow.  We  diminished  the  horizontal  response,  at  the  cost  of  a 
quicker  descent,  by  reducing  the  width  of  the  tail  array  (there  were  no  significant 
changes  in  vehicle  tilt  or  vibration).  Our  desire  for  a  fairly  steady  fall  rate  precludes 
lessening  the  negative  buoyancy  to  compensate  for  the  weaker  drag.  Weighing  the 
important  benefits  to  the  ACM  measurements  and  processing  against  the  minor  loss 
in  resolution  of  temperature  microstructure,  we  conclude  that  the  small-tail  configu¬ 
ration  is  a  valuable  advance  in  MSP’s  design. 
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Appendix  A:  SENSORS  AND  DATA  PROCESSING 


1.  ACM:  Noise  and  Aliasing 

The  acoustic  current  meter  detects  the  relative  velocity  of  water  moving  between 
its  paired  transducers.  At  intervals  of  250  /is,  each  pair  simultaneously  transmits 
2.75-MHz  pulses  of  75- pis  duration,  whose  phases  upon  arrival  at  opposing  sides  differ 
owing  to  the  water  motion.  Received  signals  are  heterodyned  with  a  frequency  156  Hz 
greater  than  that  of  the  transmission,  and  the  output  is  bandpassed  to  retain  only 
the  156  Hz  difference  component.  Further  signal  processing  yields,  at  about  312  Hz, 
voltages  proportional  to  phase  differences  and,  thereby,  to  velocities  (Brown  and 
Lawson,  1980).  This  ACM  output  is  low-passed  through  a  10-Hz  Chebyschev  filter, 
subsamplcd  at  62.5  Hz,  and  digitized  with  16  bits.  The  static  sensitivity  of  the  ACM 
is  kc2/df ,  where  k  =  0.025  is  a  circuit-related  calibration  factor,  c  is  the  speed  of 
sound  (in  meters  per  second),  d  is  the  transducer  separation  (0.2  m),  and  /  is  the 
carrier  frequency  (2.75  MHz);  nominally,  using  c  —  1500  m  s_1,  the  sensitivity  is 
0.1023  (m  s_1)/V.  We  used  this  nominal  value  until  the  Florida  Straits  cruise  before 
which  we  calibrated  the  ACM  by  repeatedly  towing  it  in  a  salt-water  tank. 

When  formed  over  short  intervals,  spectra  from  either  ACM  axis  hit  a  noise 
floor  near  10“8  (m  s-1)2  Ilz"1  and  drop  sharply  at  10  Hz  owing  to  the  low-pass  filter 
(Figure  Al).  Narrow-banded  oscillations  of  drifting  frequency  show  up  as  spectral 
spikes,  such  as  those  at  1.6  and  9  Hz  in  Figure  Al.  They  rise  around  an  order  of 
magnitude  above  the  noise  floor,  although  their  amplitudes  differ  between  axes  and 
from  cruise  to  cruise.  We  reproduced  the  behavior  by  cooling  one  of  the  oscillator 
crystals  and  forming  spectra  as  it  wanned  to  room  temperature.  Circuit  modifications 
were  then  made  that  reduced  the  spiking  and  lowered  the  noise  floor.  In  so  doing,  the 
carrier  frequency  (/)  was  slightly  lowered,  necessitating  the  empirical  recalibration. 
Because  of  their  drifting  frequency,  the  spikes  are  smeared  out  in  ensemble  spectra 
formed  over  large  intervals  and  multiple  drops,  effectively  doubling  the  level  of  the 
noise  floor. 

The  effects  of  digitally  filtering  and  subsampling  the  ACM  data  are  shown  in 
Figure  A2,  which  tracks  velocity  spectra  through  four  stages  of  processing. 

1,  The  initial  spectrum  is  formed  from  linearly  calibrated,  fully  sampled  ACM 
data;  then  the  average  fall  rate  over  the  indicated  drops  and  pressure  range  is 
used  to  convert  from  frequency  to  wavenumber.  For  each  ACM  axis,  spectra 
computed  for  0.08-MPa  (8-m)  segments  spanning  the  pressure  range  are  ensem¬ 
ble  averaged  to  form  the  autospectrum.  Displayed  in  Figure  A2  is  the  total 
spectrum,  the  sum  of  the  two  component  spectra. 
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Figure  Al.  Drop  0095,  single-axis  spectra  from  accelerometer  (top)  and  ACM  (bot¬ 
tom).  Each  spectrum  is  the  ensemble  average  of  successive  10-s  segments  (about  2.5  m) 
spanning  the  60-s  interval  from  6.8  to  6.95  MPa. 


2.  Smoothing  the  data  with  a  2-kPa  (0.2-m)  Bartlett  window  attenuates  the  spec¬ 
trum  past  1  cprn  and  halves  therms  noise.  The  effect  is  equivalent  to  multiplying 
the  un filtered  spectrum  by 


//bw  (*)  = 


(  sin(27r/rZ,/4)\  4 
\  27 vkL/4  ) 


where  /,  =  0.2  in  is  the  full  width  of  the  window  and  k  is  in  cycles  per  meter. 
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Figure  A2.  Example  of  effects  of  filtering  and  subsampling  on  ACM  velocity  spec¬ 
tra,  five  Tropic  Heat  2  drops  at  11.5°N.  Autospectra  of  the  two  orthogonal  components 
are  summed  to  yield  the  displayed  spectra  of  total  horizontal  variance.  Spectra  rep¬ 
resent  (1)  linearly  calibrated  62.5-Hz  data;  (2)  data  from  (1)  filtered  with  a  20-kPa 
Bartlett  window;  (3)  ocean  velocity  data  on  a  10-kPa  grid  (nwat);  (4)  spectrum  from  (3) 
compensated  for  the  20-kPa  Bartlett  filter. 


3.  Subsampling  the  filtered  data  onto  a  1-kPa  (0.1-m)  grid  folds  variance  past 
the  Nyquist  wavenumber  of  5  cpm  to  lesser  wavenumbers,  leading  to  noticeable 
aliasing  from  3.5  to  5  cpin.  This  spectrum  represents  the  1-kPa  ocean  veloci¬ 
ties  (t’wat)  obtained  by  removing  vehicle  motion  from  geoinagnetically  resolved 
ACM  data.  The  3.5-MPa  pressure  range  is  divided  into  1-MPa  (100-in)  seg¬ 
ments,  with  successive  segments  overlapped  by  50%.  For  each  segment.,  the 
velocity  data  (1000  points)  are  first-differenced,  demeaned,  and  llann  filtered. 
Singleton's  (1969)  method  is  used  to  take  the  Fourier  transform,  and  the  result¬ 
ing  autospectra  are  corrected  for  the  first-differencing  and  for  the  variance  lost 
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from  the  Hanning  window.  The  total  velocity  spectrum  is  the  sum  of  the  east 
and  north  autospectra.  Finally,  spectra  from  all  1-MPa  segments  are  ensemble 
averaged. 

4.  The  last  step  divides  the  previous  spectrum  by  IJ^Vi,(k)  to  recover  the  variance 
filtered  out  by  the  Bartlett  window. 

The  initial  (1)  and  final  (4)  spectra  coincide  from  1.5  to  3.5  cpm.  The  upward  hooking 
of  the  final  spectrum  at  its  highest  wavenumbers  emphasizes  its  contamination  by 
aliasing.  They  differ  at  lower  wavenumbers  because  of  the  conversion  from  relative  to 
ocean  velocity  and  also  because  of  the  difference  in  windowing  and  associated  Hanning 
weights:  segment  sizes  for  raw  velocity  were  0.08  MPa  while  those  for  ocean  velocity 
were  1.0  MPa. 

Spectra  exhibited  earlier  in  this  document  were  formed  using  the  processing  steps 
described  for  stages  (3)  and  (4);  shear  spectra  were  obtained  by  multiplying  velocity 
spectra  by  the  square  of  the  radian  wavenumber,  (27 rk)2. 

2.  Accelerometers  and  Q  rantization  Noise 

Prior  to  1989,  Sundstrand  Q-flex  model  QA-1100  accelerometers  were  used  in 
MSP  with  nominal  calibrations.  They  were  replaced  by  model  Q A- 1400  sensors, 
each  with  its  own  manufacturer-supplied  calibration.  There  is  little  difference  in 
the  performance  of  the  two  models.  Quantization  sets  the  noise  floor  at  around 
5  x  1 O'9  (in  s-2)2  Hz'1  (Figure  Al).  The  rins  noise  of  2.6  x  lO'11  m  s'2  (Table  3) 
results  from  muliplying  this  spectral  floor  by  the  Nyquist  frequency  of  12,5  Hz.  This 
quantization  noise  exhibited  by  our  16-bit  A/D  is  thrice  the  level  of  [q1  j\2)— where 
q  is  the  resolution — predicted  by  Bendal  and  Piersol  (1971).  (Note  that  the  125-Ilz 
temperature  and  conductivity  channels  are  also  quantization-limited.)  To  avoid  con 
tamiuating  ACM  velocities,  noise  is  filtered  out  of  the  small-scale  (k  >  1  cpm)  ac¬ 
celerometer  data  before  tilt  and  oscillation  corrections  are  computed. 

3.  Sensor-to-Instrument  Transformation 

An  orthogonal  coordinate  system  is  defined  on  the  vehicle  as  a  reference  for  the 
sensor  systems.  This  right-handed  ijk  frame  has  the  vertical  k  axis  up  the  geometric 
centerline,  the  i  axis  radially  outward  toward  an  inked  line  on  the  main  tube,  and 
the  j  axis  90°  anticlockwise  from  the  i  axis  (looking  from  above).  The  horizontal 
components  of  each  sensor  system  define  ab  coordinates,  with  the  magnetometer  also 
measuring  along  a  vertical  c  axis.  The  magnetometer  and  accelerometer  axes  are 
closely  aligned  with  the  instrument  frame,  while  the  ACM  and  ECM  are  rotated  in 
the  ij  plane.  Zero  offsets  (bias)  and  uncertainties  in  sensitivity  and  alignment  of  the 
sensor  axes  are  considered  as  data  are  mapped  to  ijk  coordinates. 
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The  drifting  offsets  in  the  ECM  data  are  estimated  with  a  low-frequency  spline, 
as  illustrated  in  Figure  5.  For  the  ACM  and  accelerometers,  zero  offsets  are  estimated 
for  each  drop  by  averaging  the  data  over  a  deep,  thick  interval  encompassing  an 
integral  number  of  vehicle  rotations.  These  averages  ought  to  be  zero,  assuming  that 
small-scale  oceanic  shear  is  reasonably  stationary  with  depth  (Hayes  et  al.,  1984). 
Therefore,  deviations  from  zero  are  attributed  to  sensor  bias. 

Treatment  of  magnetometer  data  relies  on  the  constancy  of  the  geomagetic  field, 
B  =  Bh  y  +  Bz  z,  during  a  drop  (although  z  is  downward  in  the  discussion  of  the 
point-mass  model,  it  is  upward  in  this  appendix,  so  Bz  <  0  at  northern  latitudes). 
Along  with  an  angle  to  orthogonalize  the  two  horizontal  axes,  offsets  and  scaling 
factors  are  found  such  that  the  adjusted  data  describe  most  closely,  in  the  least- 
squares  sense,  a  circle  of  radius  B y.  For  the  vertical  axis,  an  offset,  a  scaling  factor, 
and  two  angular  adjustments  (described  below)  are  found  such  that  the  adjusted  data 
deviate  minimally  from  B/. 

The  transformation  of  data  from  the  measured  horizontal  components  to  the 
instrument  frame  consists  of  three  steps.  In  the  following,  (da,  db)  refers  to  the  original 
data  for  any  of  the  sensor  systems  (ACM,  ECM,  accelerometers,  or  magnetometer). 
First,  the  measurements  are  adjusted  using 

da>  —  sa(da  -  ou)  and  dy  —  S(,(<4  -  ob) 

to  apply  zero  offsets  (oa,ob)  and  scaling  factors  (sQ,s/,).  Next,  with  the  a  axis  as  the 
reference  so  that  dan  =  d„/, 

<4"  =  [dai  sin  +  <40/  cos  6t/wu. 

orthogonalizes  the  data  for  an  angular  deviation  A V’ortii  of  the  b  axis  (for  magnetometer 
and  ECM,  only).  Finally,  given  the  angle  Vw  from  the  MSP  i  axis  to  the  sensor 
a  axis. 


d,  ~  da"  cos  Vw  -  <4"  sin  i/’imt 

dj  =  dan  sin  V’iusi  +  dh'<  cos  ip-ulsl 

completes  the  conversion  to  instrument  coordinate's  (accelerometers  and  magnetome¬ 
ter  assume  i/’illsl  —  ()). 

The  vertical  magnetometer  measurements  are  offset  and  scaled,  and  the  resulting 
are  oi  thogonalized  via 

bc»  =  [/>c,  —  cos  by\.  -|-  />(,<-  sin  di/y)  «in  /  cos  t>0c  , 

where  the  e  axis  is  at  angle  b0c  to  the  k  axis  and  its  projection  onto  the  ij  plane  is 

at  angle  <Si/y  to  the  a  axis.  The  last  step  is  trivial,  b ^  —  bc". 
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4.  Instrument-to-Geomagnetic  Transformation 

The  orientation  of  the  instrument  relative  to  the  geomagnetic  frame  can  be 
expressed  in  terms  of  rotations  through  three  Euler  angles  (Wrigley  et  ah,  1969). 
Given  reliable  magnetometer  ( bi,bj,bk )  and  accelerometer  (a,. a,, a/,.)  measurements, 
these  Euler  angles  can  be  determined  and  used  in  constructing  an  orthonormal  matrix 
that  transforms  data  from  the  ijk  frame  to  the  xyz  frame  (J.  H.  Dunlap,  APL-UW, 
“Coordinate  Transformation  for  MSP  Data,”  September  1985).  The  ECM  processing 
employs  this  technique,  one  result  of  which  is  the  removal  of  all  tilt  effects  from  the 
accelerometer  and  magnetometer  data  (i.e.,  ax  =  =  bx  =  0,  a,  =  <7,  bv  =  B^,  and 

K  =  BZ). 

For  the  ACM  processing,  geomagnetically  resolved  tilt  data  are  needed  so  that 
corrections  for  platform  oscillation  can  be  made.  Therefore,  an  alternate  approach 
must  be  taken  in  which  the  rotational  orientation  of  the  vehicle  is  determined  solely 
from  the  horizontal  magnetometer  data.  The  task  is  to  find  the  angle  <j>  from  the 
eastward  x  axis  to  the  instrument  i  axis. 

Before  4>  is  computed,  the  tilt  contributions  in  b,  and  b3  are  removed.  Otherwise, 
noise  or  uncertainty  at  the  order  of  the  tilt  angles  (less  than  1°)  is  generated  in  the  ro¬ 
tation  angles.  The  tilt  along  the  instrument  axes  is  estimated  from  the  accelerometer 
data, 

Qi  —  ai/ g  and  Q3=-a3/g, 


and  is  then  used  to  adjust  the  horizontal  magnetometer  measurements, 


b,>  =  6,  —  Bz  sin  6i  and 


bji  —  bj  —  Bz  sin  0j  . 


The  rotation  angle  is  such  that 

6,< 


sin  = 


(b*  +  bj,yi2 


and 


cos  <f>  = 


by 


(bl  +  b],yn  ■ 


Now.  the  ACM,  accelerometer,  and  uncorrected  magnetometer  data  are  transformed 
to  the  geomagnetic  frame  via 


dx  =  d,  cos  <b  —  dj  sin  <j> 
dy  =  d,  sin  d>  +  dj  cos  <p 
dz  —  dfc , 


in  which  tilt  contributions  remain. 

Depending  on  the  method — the  Euler  technique  or  the  use  of  horizontal  magne¬ 
tometer  data  with  or  without,  tilt  effects — the  estimated  rotation  angle  varies  slightly. 
Owing  to  the  small  tilts  of  MSP,  the  differences  among  the  methods  are  too  small  to 
be  significant  in  geomagnotically  resolving  the  ECM  or  ACM  data. 
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5.  ECM 

The  conversion  from  the  ECM  measurements  (cr,  c„)  in  microvolts  to  velocity  in 
meters  per  second  is 

1000 

•{),.  =  - e„ 

1  Uzdc(l  +  Cl)  v 

-1000 

Vy  ~  liz  de(\  T  C|)  Cr 

where  the  geomagnetic  field  is  in  nanoleslas,  the  electrode  separation  de  is  0.42  m, 
constants  c-i  and  Ci  (0.97  and  -0.03)  are  related  to  the  approximation  of  MSP  as 
an  electrically  insulating  prolate  spheroid  (Sanford  et  al.,  1978),  and  the  fall  rate  is 
W=Wmip>0. 

6.  Tilt  Angles  from  Magnetometer  Data 

In  the  ACM  processing,  the  neglect  of  MSP’s  small  tilt  angles  places  the  trans¬ 
formed  data  in  coordinates  that  are  slightly  tipped  from  the  g<x>magnetie  frame.  The 
remnant  accelerometer  signals  allow  correction  of  the  ACM  data  for  till  oscillation  as 
well  as  position.  Tilt-induced  signals  also  remain  in  the  magnetometer  data.  We  con 
sider  these  signals,  (br,  bv,bz),  for  small  tilt  components  0X  and  0V.  With  tilt  confined 
to  the  yz  plane, 


.MLlfil  w 

Bz(U-  c,) 


be  =  0 

bv  =  Liz  sintfv  -f  Bn  cos  0V 
b2  —  Bz  cos  0V  —  Hu  sin  0y  . 

When  tilt  in  the  xz  plane  is  added,  some  of  b,  projects  onto  the  x  axis,  with 

bx  =  Bz  cos  0V  sin  0t.  —  Bn  sin  0y  sin  0r 

by  =  Bz  sillily  -I  Bn  cos 0y 

b:  —  Bz  cos  0y  cos  0X  —  Du  sin  0U  cos  0X 

as  the  approximate  result  (the  same  result  occurs,  to  first  order  in  the  small  angles, 
when  considering  the  tilts  in  reverse  order  or  when  following  the  exact  and  more 
lengthy  analysis).  Therefore,  for  small  0X  and  0U  and  from  the  above  expressions,  we 
can  estimate 

0T  -  y  and  0V  -  — — - — —  =  — -  [radians]  (Al) 

Uz  Hz  Dti 

from  the  lilt-affected  magnetometer  data. 
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Appendix  B:  DATA  PROCESSING  FOR  AIRFOIL  PROBES 


1.  Sensitivity  and  Calibration  of  Airfoil  Probes 

An  airfoil  probe  responds  to  velocity  fluctuations  that  push  on  its  sensitive  tip, 
deflecting  the  piezoelectric  bimorph  beam  contained  within  (Osborn  and  Crawford, 
1080;  Oakcy,  1982).  The  induced  voltage  is  approximately  proportional  to  the  forcing 
normal  to  the  beam.  As  discussed  in  Appendix  C  regarding  the  vehicle  response, 
potential  (low  theory  predicts  a  linear  pressure  force  across  a  slender,  axisymmetric 
body  for  flows  in  which  the  lateral  component,  u,  is  small  relative  to  the  axial,  w.  The 
total  pressure  force,  F  oc  vw,  is  distributed  evenly  along  the  paraboloidal  tip.  We 
use  Oakey’s  method  for  calibration,  oscillating  the  probes  in  a  laminar  jet  of  speed 
xv.  Peak -to- peak  voltages  are  (it  against  Av  —  in  sin  AO,  the  range  of  lateral  velocity 
implied  by  the  known  oscillatiou  magnitude,  to  determine  the  calibration  factor,  Sv 
(in  volts  per  meter);  inferring  w  from  the  measured  pressure  head  (w2  =  2 gh )  accounts 
for  the  odd  units.  Now,  with  MSP  falling  at  speed  IT,  a  signal  of 

5vH  [V] 


E  = 


‘1<J 


should  result  from  a  turbulent  fluctuation  velocity  of  v.  With  increasing  angle  of 
attack,  however,  cross  drag  forcing  also  contributes  significantly,  and  therefore  the 
above  formulation  will  underestimate  the  induced  voltage  (Osborn  and  Crawford, 
1980). 


Airfoil  signals  are  conditioned  before  they  are  presented  to  the  A/D  converter. 
The  voltage  is  boosted  by  a  charge  amplifier  (in  proportion  to  the  probe’s  effec¬ 
tive  capacitance,  C,)  and  a  gain  adjusting  amplifier,  higher-frequency  (smaller- scale) 
fluctuations  aie  emphasized  by  .i  differentiating  circuit,  and  high-frequency  noise 
is  suppressed  by  a  flOTlz  four -pole  Chebyschev  filter.  The  probe  itself  has  depen¬ 
dencies  both  spatial,  as  the  tip  smooths  out  small-scale  fluctuations,  and  temporal, 
as  the  beam  equilibrates  to  attenuate  slow  variations.  Altogether,  the  transforma¬ 
tion  from  a  velocity  spectrum,  (n>2  -s~2  Hz-1],  into  a  raw  voltage  spectrum, 

4»iiw ’(/)  [V'J  Hz-1]  is  given  by 


U'L(f)  "m ,n(/W  A2  =  4®/)  , 

where  S  —  SvW/‘2g.  The  combined  electronic  response.  //£ L,  and  an  example  of 
a  raw  spectrum  (computed  with  the  robust  method  described  below)  are  shown 
in  Figure  Bl.  For  our  probe  response,  H^lfi(k),  we  use  the  formulation  of  Nin- 
nis  (1984),  a  quarlic  fit  to  measurements  valid  for  k  <  100  cpm;  we  do  not  con¬ 
sider  the  probe’s  temporal  response,  as  it  is  too  slow  to  affect  microscale  mea¬ 
surements.  Conversions  to  wavenumber  and  to  spectral  density  per  wavenumber 
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Figure  Bl.  Drop  0267,  raw  voltage  spectrum  for  airfoil  probe  #  1  computed  with 
the  robust  method  from  3.368-3.392  MPa  (K'  =  0.3241  m  s-1).  The  combined  transfer 
function  of  the  electronics,  // £L,  slopes  upward  to  25  Hz  owing  to  the  differentiating 
circuit  and  downward  thereafter  owiug  to  the  low  pass  filter. 


are  k  =  f/W  and  $^m\k)  —  The  shear  spectrum  is  computed  as 

*[£m)(k)  =  Altogether, 


tin 


(*)  = 


H7 

UlH(k) 


[s-2  cpm'1] 


where  H^H(k)  -  (S/2nk)2  //|L(/)  H^m(k)  and  /  =  kW.  Using  this  expression, 
the  raw  spectrum  of  Figure  Bl  is  transformed  to  the  shear  spectrum  of  Figure  B2. 
Also  shown  is  the  probe  response,  and  the  combined  response,  (at  this  fall 

rate).  The  nearly  constant  value  of  //|H  at  k  <  20  cpm  indicates  that  the  conditioned 
measurements  represent  shear  rather  than  velocity.  The  approximate  relationship  is 


m 


G0(C,/Cr)SvW* 

4t/o5 


1 

|//sh(/~/o)I 


[V] 
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Figure  B2.  Drop  0207.  airfoil  vs  Nasmyth  shear  spectrum.  Bottom  panel  shows 
shear  spectrum  for  airfoil  probe  #  1  for  the  robust  window  from  3.IMW  3. .'192  MPa.  The 
integration  cutoff  wavenumber,  A>  =  19.4  cpm,  is  indicated.  Dashed  line  is  the  Nasmyth 
spectrum  for  <  4.25  X  10~a  W  kg-1  and  //  =  1.22  x  lO"1’  up  s_1.  Uppoi  panel  shows 

the  Niunis  probe  response,  1N  (light  line),  and  the  combined  transfer  function  from 
shear  to  measured  voltage,  Jl^ j  (heavy  line). 


where  6‘0  =  9.7  is  the  circuit  gain  at  /0  =  1  11/,  ami  (',/(' j  fs  1  is  the  further,  probe- 
dependent  gain. 

2.  Spectral  Processing 

The  periodogram  processing  generates  dissipation  rates  on  a  5  k Pa  (0.5-in)  grid. 
Two  values  —one  for  each  airfoil  probe-  -are  computed  for  each  point.  Eacli  estimate 
is  based  on  data  from  a  window  of  L  =  5  kl'a  centered  about  the  grid  point,  yielding 
no  overlap  between  windows.  Since  the  raw  data  are  sampled  at  125  IIz,  the  number 
of  points  in  a  window  varies  with  fall  rate.  After  the  data  are  demeaned,  a  Han 
ning  window  (full  cosine  taper)  is  applied.  Zeroes  are  appended  to  the  record  until 
attaining  a  valid  l'FT  size  (an  even  integer  factorable  into  2,  3,  and  5).  The  power 
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spectrum  is  formed  from  the  transform  of  the  data.  After  this  raw  frequency  spec¬ 
trum  is  corrected  for  the  variance  filtered  out  by  the  time-IIanuixig,  it  is  converted  to 
a  wavenumber  spectrum  of  shear. 

The  robust  processing  estimates  dissipation  rates  on  a  20-kPa  (2-m)  grid.  It  uses 
data  windows  of  L  =  24.24  kPa  centered  about  each  grid  point.  The  slight  overlapping 
of  windows,  combined  with  the  spectral  processing  within  them,  evenly  weights  all 
airfoil  data  in  the  drop.  Multiple  spectra  are  computed  within  each  window  on 
segments  of  length  £'«L/4  =  G*06  ltl’a.  The  actual  size,  which  varies  with  fall  rate,  is 
chosen  such  that  the  number  of  raw  data  points  is  a  valid  FFT  size.  The  raw  spectrum 
for  each  segment  is  formed  using  the  same  steps  as  in  the  periodograui  processing  (but 
without  zero- padding).  Segments  are  overlapped  by  approximately  70%,  yielding  9 
to  11  spectra  per  window  (Figure  133).  As  spectral  values  at  each  frequency  are 
combined  to  form  the  ensemble  spectrum  for  the  window,  m-estimators  are  used  to 
heavily  downweight  high  outliers  (Cliave  et  ah,  1987).  This  technique,  implemented 
by  W.  R.  Hess  (personal  communication,  1989),  objectively  and  reproducibly  rejects 
contamination  from  spikes  caused  by  plankton  impacts  or  electrical  transients  in 
MSP’s  circuitry.  However,  the  algorithm  fails  when  spiking  occurs  too  often  in  a 
window.  Also,  if  the  amplitude  of  the  true  signal  changes  abruptly,  the  window’s 
spectrum  may  be  underestimated  as  the  higher  variance  segments  are  mistakenly 


time  /  s 

0  5  10  15  20 


Figure  B3.  Windowing  lor  robust  processing  of  airfoil  data,  assuming  VT  =  0.40  m  s_1. 
Windows  art1  2.424  in  wide,  centered  at  2,  4,  and  0  m.  Each  window  consists  of  10 
segments,  0.000  m  long  and  overlapped  70%.  Shape  of  segments  indicates  weighting  of 
Hanning  windows. 
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downweighted.  This  effect  is  slight  for  the  airfoil  data,  but  it  can  be  severe  in  robust 
processing  of  temperature  gradient  microstructure. 

The  two  methods  produce  similar  results  for  sections  of  clean  data,  judging 
from  dissipation  rates  averaged  over  0.1  MPa  or  more.  We  also  find  good  agreement 
when  comparing  spectra  ensemble  averaged  over  several  consecutive  windows;  larger 
windows  are  used  in  this  case  than  in  generating  dissipation  rates  to  provide  a  wider 
overlap  with  ACM  results  (such  as  in  Figure  17).  There  is  a  tendency  for  robust 
dissipation  rates  or  spectra  to  be  slightly  (up  to  10%)  lower  than  the  periodogram 
results. 

3.  Dissipation  Rates 

For  isotropic  turbulence,  c  =  7  .bv(du2/dz)2,  where  v  is  the  kinematic  viscosity 
and  u2(z)  are  the  horizontal  cross-stream  velocity  fluctuations.  We  estimate  the  shear 
variance  by  integrating  ^a^m\k)  over  an  appropriate  wavenumber  band.  The  size  of 
the  segment  chosen  in  computing  the  spectrum  sets  the  lower  limit,  k0  (as  well  as 
the  wavenumber  interval,  A k  =  k0).  Regarding  the  upper  limit,  ku ,  it  should  be 
high  enough  so  that  most  of  the  variance  has  been  captured  but  not  so  high  that 
noise  significantly  contaminates  the  integration.  The  spectral  noise  floor  is  probe- 
dependent  and  sometimes  varies  during  a  drop.  However,  the  lower  the  dissipation 
(and  corresponding  spectral  levels),  the  sooner  the  noise  floor  is  reached,  and  the 
sooner  the  integration  must  be  cut  off.  We  use  Oakcy’s  analytic  representation  of 
the  Nasmyth  turbulent  shear  spectrum  to  guide  our  selection  of  ku.  This  model 
spectrum  has  a  constant  shape,  but  shifts  in  wavenumber  and  amplitude  with  varying 
Kolmogorov  wavenumber,  k,  =  (c/V3)1/M,  and  spectral  level,  A,  =  (a/5)1^4. 

Summing  the  data  spectrum  from  ku  up  to  10  cpm  generates  an  initial  estimate, 
ciocpm-  determined  by  finding  the  model  spectrum  that  contains  the  same  variance  as 
the  data  over  this  limited  band.  The  trial  cutoff  of  ku  —  0.088  k„  (ku  in  cycles  per 
meter)  corresponds  to  the  wavenumber  before  which  resides  90%  of  the  variance  in 
the  model  spectrum  for  If  ku  <  9  cpm,  it  is  set  to  9  cpm  (this  occurs  when 

(  <  2  x  HI- 10  W  kg-1);  if  k\  >  100  cpm,  it  is  set  to  100  cpm,  the  upper  limit  of  tin; 
probe’s  resolution  (this  occurs  when  c  >  10_r>  W  kg-1).  Summing  the  data  spectrum 
from  A’0  to  ku  produces  an  improved  estimate,  <„t,  A  new  ku  is  determined  from 
tet,  and  the  spectrum  is  again  summed  to  yield  the  final  value  of  c.  No  attempt  is 
made  to  compensate  for  unresolved  variance  (i.e.,  below  1  cpm  or  past  100  cpm), 
since  dissipation  rates  encountered  by  MSP  seldom  exceed  10~5  W  kg-1.  At  each 
grid  point,  two  estimates  (one  from  each  probe)  are  generated.  If  they  are  close  (e.g., 
within  a  factor  of  4),  they  are  averaged  together;  otherwise,  the  lower  value  is  chosen 
since  noise  and  spiking  lead  to  inflated  estimates. 

In  Figure  B4a  we  compare  a  range  of  data  spectra  with  their  corresponding 
Nasmyth  spectra.  Each  data  spectrum  is  formed  by  ensemble  averaging  0.01-MPa 
periodogram  spectra  (50%  overlapped)  over  intervals  of  0.1  to  0.2  MPa,  within  which 
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Figure  B4.  Comparison  of  ensemble-averaged  airfoil  spectra  (heavy  lines)  with  (a) 
corresponding  Nasmyth  spectra  (light  lines)  and  (b)  corresponding  single-axis  ACM 
spectra  (dashed  lines).  Drop  numbers,  pressure  range,  t,  and  T  for  the  spectra  are, 
from  bottom  to  top: 


034G,  4.03- 

4.17 

MPa, 

3.2 

X 

10- 

-n  W  kg"1, 

10.3°C, 

0387,  4.65- 

-4.70 

MPa, 

1.6 

X 

io- 

'10  W  kg-1, 

8.8°C, 

0348,  9.35- 

-9.47 

MPa, 

4.2 

X 

io- 

'9  W  kg-1, 

4.7°C, 

0356,  2.79- 

2.91 

MPa, 

1.2 

X 

io- 

-7  W  kg-1, 

11.9°C, 

0267,  4.83- 

-4.93 

MPa, 

1.6 

X 

io- 

-6  W  kg-1, 

7.8°C. 

Arrows  in  (a)  indicate  average  cutoff  wavenumbers  for  t  integrations. 
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there  is  minimal  variation  of  dissipation  rates.  Also,  spectra  from  the  two  probes 
are  averaged  together,  since  they  agreed  well  in  all  cases.  The  Nasmyth  spectrum 
is  formed  using  the  mean  of  the  0.02-MPa  robust  dissipation  estimates  spanning  the 
interval.  All  but  the  highest  spectrum  hit  the  noise  floor  by  100  cpm,  with  the  lowest 
spectrum  barely  above  the  noise  at  k  <  5  cpm.  Integration  of  such  low  spectra  (out 
to  9  cpm)  establishes  our  dissipation  noise  level  of  3  to  10xl0-11  W  kg-1.  At  low 
dissipation  rates,  the  variance  of  the  inertial  subrange  is  concentrated  below  1  cpm, 
beyond  the  limit  for  the  airfoil  probes;  oftentimes,  there  is  no  inertial  subrange, 
with  motions  becoming  anisotropic  owing  to  buoyancy  effects  (e.g.,  for  t/uN2  <  25). 
As  rates  rise  above  2  x  1G-1U  W  kg-1,  turbulent  shear  spectra  are  increasingly  well 
resolved.  The  measured  and  model  spectra  agree  decently,  with  minor  differences  at 
low  wavenumber  and  in  the  location  of  the  viscous  rolloff. 

Combining  uncertainties  of  fall  rate,  probe  calibration,  temperature  sensitivity, 
and  spectral  processing  and  integration,  we  estimate  that  our  measured  dissipation 
rates  are  within  a  factor  of  2  the  actual  values.  However,  comparisons  with  the 
ACM  suggest  additional  uncertainty  of  an  unknown  nature  that  might  overwhelm 
the  above-mentioned  sources. 


4.  Airfoil  Probes  vs  ACM 

A  unique  feature  of  MSP  is  that  it  measures  velocity  micros tructu re  with  two 
separate  devices.  Tin*  ACM  averages  motions  across  20  cm  paths,  positioned  2  cm 
ahead  of  the  point-like  measurements  of  the  airfoil  probes.  The  alignment  of  the 
airfoils  differs  from  that  of  one  ACM  axis  by  only  20  25°.  We  compare  shear  spectra 
from  this  axis,  uncorrected  for  platform  motion,  with  from  the  airfoils.  The 

spectra  often  agree  well  from  1  to  10  cpm  (or  to  where  the  ACM  becomes  noisy); 
coherence  of  the  data  can  be  even  more  broadbanded,  from  0.1  to  20  cpm  (at  most). 


With  such  comparisons,  we  demonstrate  the  general  success  of  the  airfoil  calibra 
tions  (when  l  CM  and  airfoil  spectra  coincide),  diagnose  contamination  by  elevated 
low-frequency  noise,  and  note  an  apparent  increase  in  airfoil  sensitivity  with  tempera¬ 
ture  and  dissipation  rate.  We  infer  contamination  when  the  spectrum  from  the  failing 
probe  is  well  above  that  of  the  ACM  at  1  cpm  and  steadily  falls  to  join  it  at  some 
greater  wavenumber  (while  the  other  probe  tracks  the  ACM  throughout);  in  more 
turbulent  sections  of  the  drop,  the  spectral  mismatch  disappears  as  the  real  signal 
exceeds  the  low-frequency  noise.  In  contrast,  variations  in  sensitivity  are  suggested 
when  the  airfoil  probes  mutually  agree  but  change  from  nearly  matching  to  greatly  ex¬ 
ceeding  the  ACM  as  conditions  vary  during  a  drop.  Most  often,  the  ACM  spectrum 
runs  nearly  parallel  to  the  higher  airfoil  spectrum,  but  sometimes  their  separation 
increases  with  wavenumber  (as  if  the  ACM  spectrum  is  attenuated). 


In  Figure  1Mb  we  plot  airfoil  spectra  against  their  ACM  counterparts.  For  the 
lowest  pair,  the  ACM  is  barely  above  its  noise  floor.  The  next  highest  show's  good 
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agreement  from  1  to  4  cpm.  Moving  up,  the  next  two  are  well  matched  in  shape,  but 
the  airfoil  spectra  are  higher  by  factors  of  1.3  and  2,  respectively.  The  highest  spectra 
differ  in  shape,  diverging  toward  higher  wavenumbers;  however,  the  ACM  spectrum 
for  tlie  next  deeper  10-m  interval  of  this  drop  parallels,  at  a  factor  of  2.2  beneath,  its 
airfoil  spectrum  (which  is  close  in  level  to  that  displayed). 

We  characterize  each  comparison  by  summing  both  spectra  from  2  to  7  cpm  and 
then  forming  the  airfoil-to-ACM  ratio.  For  c  <  10"9  W  kg-i,  we  stop  at  5  cpm  to 
avoid  ACM  noise  (and  at  4  cpm  for  e  <  3  x  lO'10  W  kg-1).  We  summarize  results 
from  the  COA11E3  cruise  in  Figure  B5,  plotting  ratios  against  temperature,  T\  and 
dissipation  rate,  c.  There  is  fair  agreement  at  lower  T  and  £,  with  ratios  varying  from 
0.5  to  2  and  occasionally  3.  Ratios  are  generally  greater  at  the  highest  T  and  c,  often 
rising  above  4  and  to  as  much  as  8.  Correlations  between  c,  I\  and  N  (or  dT/dz) 
confound  attempts  to  distinguish  the  source  of  the  observed  behavior.  The  ratio  does 
rise  with  t  for  T  «  29°  and  with  T  for  c  «  3  x  ll)~9  W  leg-1.  Comparisons  from  other 
MSP  cruises  show  patterns  consistent  with  those  of  Figure  135. 

Most  of  our  published  shear  spectra  are  computed  within  the  thermocline,  where 
T  <  20°C  and  c  rarely  exceeds  10~7  W  kg-1;  thus,  ensemble-averaged  ACM  and  airfoil 
levels  appear  well  matched  between  1  and  3  cpm.  It,  was  not  until  we  focused  on  the 
warm,  turbulent  waters  atop  tin*  tropical  Pacific  thermocline  that  we  noticed  the 
spectral  mismatches. 

Variations  in  airfoil  or  ACM  sensitivities  can  lead  to  some  spectral  separation. 
When  an  airfoil  probe’s  sensitivity  (or  effective  capacitance,  C,)  is  greater  than  sup¬ 
posed,  use  of  the  lower,  supposed  value  will  boost  the  computed  shear  spectrum.  In 
contrast,  if  the  ACM  sensitivity  is  greater  than  supposed,  the  computed  spectrum  will 
be  too  low.  Indeed,  airfoil  sensitivity  can  increase*  with  temperature  and  turbulent 
velocity  (Osborn  and  Crawford,  1980)  and  ACM  sensitivity  does  go  up  with  sound 
speed  (and,  therefore,  with  temperature  for  MSB’s  depth  range).  However,  these 
effects  account  for,  at  most,  a  factor  of  2  difference  between  airfoil  and  ACM  spectra. 
Furthermore,  the  observed  behavior  lacks  consistency;  for  instance,  spectral  ratios 
at  a  given  temperature  and  dissipation  rale  can  vary  considerably  from  case  to  case 
(even  in  a  single  drop).  Other  possible  contributors  include  temperature  sensitivity 
in  the  circuits  of  the  airfoils  (probably  weak)  and  the  ACM  (complex  and  diflicult 
to  quantify)  and  attenuation  across  the  ACM  path  (difficult  to  characterize).  One 
point  in  support  of  the  contention  that  the  airfoil  is  responsible:  when  their  spectral 
levels  are  reduced  to  agree  with  those  of  the  ACM  (in  Figure  B4  for  instance),  the 
lowered  airfoil  spectra  and  Nasmyth  spectra,  recomputed  using  the  correspondingly 
lower  dissipation  rates,  tend  to  exhibit  a  better  match  than  the  originals,  particularly 
in  the  high  wavenumber  rolloff.  If  the  airfoils  are  the  cause,  we  are  still  uncertain 
whether  the  behavior  is  confined  to  MSP  probes  and  circuits,  to  all  probes  of  our 
design  and  manufacture,  or  to  all  piezoelectric  airfoil  probes. 
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Figure  B5.  Ratios  of  airfoil- to- ACM  shear  variance  from  C0AIIE3  drops.  Shear 
spectra,  formed  for  0.1  to  0.3  MPa  intervals,  are  integrated  from  2  to  7  cpm  to  compare 
variance  (to  av  id  ACM  noise,  integrations  end  at  4  or  5  cpm  for  low-level  spectra). 
Ratios  are  indicated  by  symbols,  which  arc  plotted  against  the  average  <  and  T  for 
each  case.  The  intervals  chosen  for  comparison  include  most  of  the  high  dissipation 
occurrences  (<  >  10-8  W  kg-1)  but  only  a  representative  sampling  of  lower  dissipation 
cases.  The  isolated  high  ratio  at  (3  x  10-9  W  kg-1  and  13°C)  is  due  to  spiky  data  in 
both  airfoil  channels. 
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Appendix  C:  VEHICLE  DYNAMICS 


1.  Horizontal  Dynamics 

We  restrict  the  horizontal  dynamics  by  assuming  MSP  remains  upright  at  all 
times.  Therefore,  6(t)  =  d9(t)/dt  =  0,  and  ali  points  on  the  vehicle  move  at  vmsp(<), 
the  velocity  of  P  from  (3).  The  dynamics  are  like  those  of  a  point  mass,  since  the 
vehicle  must  accelerate  uniformly  regardless  of  where  it  is  forced.  The  water  trapped 
under  the  skin  moves  along  with  the  solid  structure,  so  we  include  it  in  the  defini¬ 
tion  of  A/,  the  total  vehicle  mass.  We  examine  the  forces  acting  on  MSP  when  the 
ACM  is  at  reference  depth  z  =  z(t),  as  defined  by  (6).  At  some  distance  s  >  0 
above  the  ACM,  the  water  encounters  MSP  with  a  (free  stream)  relative  velocity 
vWat(T  —  s)  —  v^p (f),  the  difference  between  (1)  and  (3).  While  the  transverse  (hor¬ 
izontal)  component 

^wat(^  •$)  ^iasp(0 

can  vary  along  the  body  (Figure  9),  we  prohibit  variation  of  the  axial  (vertical) 
component 

Ww* (z  -  s)  -  vwo 

by  assuming  that  VFwat  is  constant  over  MSP’s  length. 

The  horizontal  pressure  forces  employed  by  the  model  are  based  on  potential 
flow  theory.  We  summarize  this  result,  following  Hayes  et  al.  (1984).  Consider  a  long, 
slender  solid  body  of  revolution  about  some  s  axis,  fixed  in  place  and  pointed  into 
a  uniform,  inviscid  flow  with  axial  and  transverse  components  w  and  v,  respectively 
(Figure  Cl).  For  a  small  angle  of  attack  a  =  tan ~1(v/w),  the  body  experiences  a 
force,  /(s),  per  unit  length  in  the  transverse  direction  of  magnitude. 

dA  dA 

f{s)  =  q0—  sin 2ct  =  pvw—  [N  m"1]  ,  (Cl) 

as  as 

where  =  p(v2  +  w2)/2  is  the  free-stream  dynamic  pressure,  p  is  the  fluid  density, 
and  A{s)  =  7 r(r(s))2  is  the  cross-sectional  area  of  the  body  (Allen  and  Perkins,  1952). 

To  apply  (Cl)  to  MSP,  we  fix  the  origin  of  the  s-axis  at  z  —  z(t)  with  s  >  0 
upwards  (Figure  9),  and  sum  f(s)  up  along  the  instrument.  Since  f(s)  —  0  where 
dA/ds  —  dr/ds  —  0,  the  constant -diameter  main  tube  should  feel  no  linear  pressure 
force.  Integration  of  (Cl)  from  the  ACM  (r  =  0,  s  —  0)  up  onto  the  front  of  the  main 
tube  {r  —  rA,,  s  as  1  m,  dr/ds  —  0)  yields  the  nose  force,  Fn\  integration  from  the 
back  of  the  main  tube  (r  rN,  s  «  3  m)  up  to  the  maximal  radius  at  the  brushes 
(r  =  rT ,  s  ~  3.7  m)  yields  the  tail  force,  Fj.  While  we  do  not  expect  a  uniform  flow 
in  either  region,  we  do  assume  that  the  effective  horizontal  flow  felt  near  the  nose 
or  tail  can  be  represented  by  a  single  relative  velocity,  Vn  or  Vj.  Specifically, 
is  some  suitable  average  of  the  horizontal  relative  velocities  L'wat(^  —  s)  —  within 
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Figure  Cl.  Pressure  force  per  unit  length  on  a  slender  body  of  revolution  whose 
axis  of  symmetry  is  at  angle  a  to  an  oncoming  uniform  flow,  v.  Cross-sectional  area  is 
vl(s)  =  7r(r(.s))2. 

the  integration  bounds.  For  the  oncoming  axial  flow',  | |  =  H7wal  >  0 

has  the  proper  sense.  Integration  of  (Cl),  under  our  stated  assumptions,  yields  the 
expressions  in  (11)  for  Fyv  and  Fj.  In  the  subsection  on  model  formulation,  we  explain 
how'  tyv  and  Vj  are  evaluated. 

Including  the  added  mass  effect,  as  described  in  the  main  text,  the  balance  of 
horizontal  forces  is 

M^  =  FN  +  FT-pVe^, 

at  at 

which  simplifies  to  (12)  for  acceleration.  Not  included  in  this  balance  is  the  transverse 
quadratic  drag 

Fdx  —  2  pCoxArv2  ,  (C2) 

where  Cdx  is  the  drag  coefficient,  v  is  the  relative  horizontal  velocity,  and  Ax  is  the 
body’s  area  normal  to  v.  In  most  instances,  drag  force  on  MSP  is  negligible  compared 
to  Fp/  and  Ft,  and  so  we  exclude  it  from  the  model.  However,  drag  on  the  main  tube 
may  be  significant  when  encountering  high  shear. 

2.  Vertical  Dynamics 

Although  computation  of  vertical  forces  is  not.  part  of  the  model,  it  is  useful  to 
express  the  associated  dynamics.  With  Ms  and  T4  as  the  mass  and  volume  of  the 
solid  structure, 
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•  Fa 

-  Fb  ~ 

Fp, 

gK 

~  9  />!'»  - 

■  \pCD,A,Wll 

(C3) 

gravity 

buoyancy 

drag 

where  Cp.  ssl  is  the  vertical  drag  eoeliicient,  aud  A ,  =  itr\-  is  the  froutal  area.  Even 
when  drag  balances  net  buoyancy,  I'q  —  1'b >  to  keep  |W/,.ei|  =  H7ltut,,  —  U7w»t  constant, 
MSP’s  absolute  fall  rate,  M7nui, ,  will  vary  as  the  water  motion,  H7will,  fluctuates  with 
depth. 


3.  Tilt  Dynamics 

We  assume  MSP  acts  like  a  forced,  damped  pendulum  so  that 


Iv~i)2  +  Hu  *h  A  i;0  =  Lp  FN  —  Lr  f'r  IN  m] 


(C4) 


where  Lp  and  Lj  are,  respectively,  the  distance  between  Fp  and  P  and  Ft  and  P. 
Coefficients  / p ,  lip,  aud  A’^,  as  described  and  roughly  estimated  in  Table  2,  are 
“effective"  in  that  they  include  the  following  contributions  arising  from  interactions 
between  the  falling  MSP  and  the  surrounding  water:  Ip  includes  added- mass  effects  as 
well  as  the  mass  distribution  of  the  vehicle  and  trapped  water;  A/.;  includes  damping 
from  oscillation- induced  pressure  forces  at  the  tail;  A#  includes  righting  from  vertical 
drag  at  the  tail  [Fp.)  as  well  as  the  separation  of  Fa  and  Fp.  The  forcing  terms  are 
of  opposite  sign  since  Fp  >  U  gives  an  anticlockwise  acceleration  while  Fj  >  0  gives 
a  clockwise  push. 


We  simplify  ((-4)  by  dividing  through  by  fp  and  by  referring  to  (11)  for  Fp  and 
Ft,  yielding 

~dP  +  ~dt.  ^  =  ^'N  ~  ^  T  tra^  s_ii  ’  (f  ’'r)) 

where  ft  -  lip/ IE  ,  ft,  =  Kp/Ip, 


Q\  =  ~T~1vr7pP  I l'l'7rd|  and  Qi  =  -r'p)p\Wtc)\  [in'1  s'1]  (CC) 

JL  ip 


and  with  Vp  and  Vj  as  specified  in  (17).  W'e  do  not  explicitly  use  (C5)  in  generating 
velocity  profiles.  Rather,  it  enters  into  the  transfer  functions,  derived  in  Appendix  E. 
with  which  we  analyze  effects  of  tilt  dynamics  on  the  ACM  and  accelerometer  mea¬ 
surements.  W'hile  (C5)  competently  models  some  aspects  of  the  observed  behavior, 
it  is  an  inadequate  substitute  for  the  tilt  as  inferred  in  (13)  from  the  accelerometer 
data. 
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Appendix  D:  MODEL  NUMERICAL  INTEGRATION 


With  i^msp (z)  constant  over  the  integral  bounds  in  (17),  we  can  express  (15)  as 

/•Zn+0.5Dh 

/  Uwal(z-s)  ds  (Dl) 

J  Ztf  —0.5  Dp] 

Hj  cZt+0.5Dt 
—  vwgii(z-s)  ds. 

Dr  J Zt— o.5Z)t 


dt-'msp  #  A  {Zn+OSDh 

j  d*  ( 77 /v  -f"  Ht  )  ^msp  ( ^ )  —  n  /  ^ 

d*.  JJ a  JZpi—O.hDpi 


+ 


The  desired  results  are  profiles  of  uinsp  and  uwat  at  evenly  spaced  depths  z_,  =  z0+j Az. 
We  define  a  function  U}  that  satisfies 


dUj 

dz 


+  (Hn  +  Ht)Uj  =  Bj 


(D2) 


in  the  interval  [zj_i,Zj],  where 

Hn 


Bi 


Dk 

Hr 


[  I(zj- 1/2  -  {ZN  +  0.5 Dn))  -  /(z;_1/2  -  (ZN  -  O.GD*))  ] 


+  ^  (  /(^--i/2  -  (Zt  +  0.5 Dt))  -  l(zj. i/2  -  (Zr  -  0.5DT))  ] 


is  the  right  hand  side  of  (Dl)  evaluated  at  the  interval  midpoint  z  =  zrJ_1/2  =  zj  —  ~ 
using  integrals 

J(Z)  =  /  tw(s) 

J  ?o 

linearly  interpolated  from  the  running  sums 


ds 


r./  ,  _  v—'  T>wa(.(z,)  +  T’watf2!- 1)  A 

M*»)  =  (L - ^ - Az . 

«= i  ^ 


(Note  that  /(Z  <zo)  =  0  and  that  6(zo)  =  0.)  With  the  initial  condition  t/j(z;_i)  — 
l,n1sp(^'-i),  this  differential  equation  for  Uj  approximates  (Dl)  for  vnvsp  over  the  inter¬ 
val,  and  thus  nIusp(z;)  ps  Uj(Zj).  The  solution  of  (D2), 


U}(z)  = 


^msp 


,zj-i ) 


Bj 


Hn  +  Ht  J 


exp[-(/fN  +  //r)(2  -  f,-i)]  + 


//a-  +  Hr 


is  evaluated  for  z  =  z}  to  update  the  vehicle  velocity,  nmsp(zJ)  =  Uj(zj).  Finally,  as 
in  (18),  we  set  nwat(zj)  —  ^msp^,?)  “b  ^rcl(^j)- 
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Appendix  E:  TRANSFER  FUNCTIONS  AND  VEHICLE  TILT 


In  the  first  section  of  this  appendix,  transfer  functions  are  derived  for  the  hori¬ 
zontal  (point-mass)  and  tilt  motion  of  MSP  as  it  responds  to  oceanic  shear.  They  are 
used  in  the  main  text  to  find  values  for  the  point-mass  parameters  and  later  in  this 
appendix  to  find  values  for  the  tilt  response  parameters.  While  not  entirely  realistic, 
the  modeled  tilt  adequately  demonstrates  the  dominance  of  forcing  on  the  response 
and  the  appropriateness  of  inferring  tilt  from  the  accelerometer  data.  The  effects 
of  tilt  position  and  oscillation  on  the  ACM  measurements  are  evaluated,  with  the 
conclusion  that  uncertainties  in  the  computed  tilt  or  in  the  correction  parameters  Hj 
and  Lv  have  little  effect  on  the  resulting  ocean  shear  spectra. 


1.  Response  Functions 

We  generated  several  of  the  parameters  listed  in  Table  6  while  discussing  vehicle 
dynamics,  interpreting  the  measurements,  and  specifying  the  vehicle  response  model. 
Their  various  effects  are  examined  most  easily  through  the  analytic  expressions  of  the 
vehicle  response  functions  in  wavenumber  space. 

We  define  the  Fourier  transform  for  water  velocity  as 

(k)=  /  I’wat  .(z)c~tk*dz 
J -oo 

with  radian  wavenumber  k  and  with  2  positive  downward.  As  indicated  in  Table  f), 
transfer  functions  A,  13,  Cr ,  Cm,  and  D  predict  how  MSP  and  its  measurements 
respond  to  the  ocean  velocity  field. 

First  we  derive  A(A:),  the  response  function  for  vehicle  velocity.  Starting  with 
the  Fourier  transform  of  (Dl), 


(?A-  +  //A-  +  Hr)  f'map W 


+ 


U  y  rZfij+Q.SDtf 

l) /V  J 7't-i  — O.SDyv 
II j  y^r+0.5Dr 

D'J'  J Z*} — 0.5i^T 


C  lk‘  Uwm(A’)  ds 
c~lkr>  Kai.[k)  ds  , 


we  divide  through  by  fwat (k).  evaluate  the  integrals,  and  simplify  using  Euler’s  in- 
dentity,  sin  a-  =  (c,r  —  c",r)/2/,  to  form 


A(k)  = 


Hn  c~'kZN 

//,v  -f  Hr  T  ik 


sine  {kDv/2)  + 


Hr  c~'kZr 
H  k  T  Hr  +  1  k 


sine  (ArDj/2)  , 


where  sine  .r  =  x  1  sin  x  . 


(El) 
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We  proceed  in  a  similar  fashion  to  derive  B(k),  the  response  function  for  vehicle 
tilt.  We  assume  a  constant  fall  rate,  Wo ,  so  that  dOjdt  —  W0dO/dz ,  and  simplify  the 
transform  of  (C5)  to  produce 

n„A.  QNe-ik*«8mc{kDN/2)-Q7  c~ikZr  sh\c(kDT/2)  —  {Qn-Qt)A(Ic)  /TViX 

B {k} - k+ikw„p,-k>w3 - •  (E2) 


We  generate  the  remaining  transfer  functions 

C,  (k)  =  1  -  A(k) 

Cm{k)  =  Cr(k)  +  \V0{B,  -ikLv)B(k) 

D(k)  =  ikW0  A(k)  +  (<j-  I,ak2W*  -  ikWg)  B(k) 


(E3) 


by  dividing  i'wa(  into  the  Fourier  transforms  of  (7),  (14),  and  (10),  for  which  we 
assume  W0  |t'mS|,0|.  Expansion  of  A(k)  in  (E3)  leads  to  (19),  with  which  we  tune 
the  point-mass  response  model. 

2.  Tuning  the  Tilt  Response 

To  complete  our  model  of  the  tilt  response,  we  must  select  values  for  parameters 
ft,  An  Qn,  and  Qj-.  We  begin  by  assuming  the  optimal  values,  given  in  Table  6,  for 
the  six  force  parameters  in  A{k).  We  could  tune  B(k)  directly  by  comparing  (E2) 
with  cross  spectra  computed  from  profiles  of  0  and  Instead,  we  prefer  to  work 
with  the  transfer  function  from  to  a„lcll,,  since  these  directly  measured  data,  and 
their  cross  spectrum,  arc  independent  of  model  computations  or  assumptions.  We 
define  this  transfer  function  as 


nm(k) 


.(*)  D(k) 


M')  Cm{k) 


(E4) 


and  note  that,  in  the  model,  it  includes  parameters  Bj,  Lv,  and  La.  We  will  also  refer 
to  the  transfer  function  from  relative  velocity  to  measured  acceleration, 


r  /  j.\ 

M*) 


^meos  (k)  D(k) 

v„i{k)  C\(k)' 


(E5) 


as  we  examine  the  effects  of  tilt  and  oscillation  corrections. 


We  generate  trial  versions  of  Hm  by  varying  /ft ,  f}0,  and  Qr  while  keeping  Bj , 
Z/„,  and  La  fixed  at  their  nominal  values  in  Table  6.  Once  we  attain  decent  agreement 
with  the  measured  cross  spectrum,  we  can  adjust  Bj  and  Lv  to  improve  the  match. 
This  tuning  process  yields  the  optimal  parameter  values  listed  in  Table  6.  In  Fig¬ 
ure  El,  we  compare  the  theoretical  forms  of  Hm  and  HT  with  the  corresponding  cross 
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Figure  El.  Comparison  of  model  transfer  functions  Hm  and  //r  with  corresponding 
cross  spectra  computed  from  PATCH  EX  da  ta.  Top  panel  is  data  coherence.  Middle  and 
bottom  panels  arc  amplitude  and  phase.  Heavy  lines  denote  the  model  and  dash-dot 
lines  denote  the  data  (95%  confidence  limits  are  shaded).  Dashed  line  in  middle  panel 
is  the  square  root  of  the  ratio  of  data  autospectra,  which  differs  in  amplitude  from  the 
data  transfer  function  where  coherence  is  small,  (a)  ACM  measured  velocity  to 
measured  acceleration;  (b)  //r,  corrected  relative  velocity  to  measured  acceleration. 


TR  9414  E3 


E82KSJS 


UNIVERSITY  OF  WASHINGTON*  APPLIED  PHYSICS  LABORATORY 


Figure  El,  coni. 


spectra  derived  from  profiles  (of  the  northward  components)  of  the  28  PATCHED 
drops.  The  top  panels  show  the  squared  coherence  of  the  data, 

1’he  match  between  the  analytic  and  measured  transfer  functions  is  beet  for 
k  <  0.3  cpm  and  degrades  at  highe.  wavenumbers.  It  is  problematic  to  find  \  set  of 
parameter  values  that  match  the  data  for  both  amplitude  and  phase,  however.  At  0.22 
and  0.42  cpm,  Hm  exhibits  large  peaks  in  amplitude  and  wide  swings  in  phase.  These 
features  correspond  to  the  notches  in  $^'IU  and  to  bands  of  minimal  coherence  in 
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the  data.  Removal  of  tilt  and  oscillation  from  the  measurements  enhances  coherence 
of  relative  velocity  and  acceleration  (i.e.,  t»re]  and  ameiul)  at  the  notch  wavenumbers. 
Comparing  //,„  and  llr  at  k  >  0.5  cpm,  the  phase  -originally  near  90° — is  disrupted 
and  coherence  is  diminished,  consistent  with  removal  of  oscillations  from  tv,,,. 

3.  Point-Mass  vs  Tilt  Response 

The  point  mass  response  predicted  by  (El),  shown  in  Figure  E2,  and  that  re¬ 
booted  by  the  cross  spectrum  of  uuup  and  profiles  are  identical  to  within  the 
accuracy  of  the  numerical  method.  The  amplitude  |/l(&)|,  which  is  essentially  the 
factor  7  discussed  with  Figure  4,  diminishes  only  after  0.33  cpm,  reflecting  MSP’s 
tendency  to  respond  at  scales  of  the  vehicle  and  larger. 


10'2  10° 

k  /  cpm 


Figure  E2.  MSP  platform  response  for  large-tail  configuration.  Light  line  is  the 
point-mass  response,  A[k).  The  remaining  curves  compare  the  combined  effects  of  tilt 
and  oscillation,  -BjWaui)0  +  L^dO/dl,  to  the  point-mass  motion,  1^,,.  Heavy  hues 
denote  the  model,  \CT(k)  —  Cm(k)]/A[k),  and  dash-dot  lines  denote  the  cross  spectrum 
of  the  PATCH  EX  data  (with  95%  confidence  limits );  dashed  hue  is  the  square  root  of 
the  ratio  of  data  aulospectra;  the  coherence  is  the  same  as  for  IIr(k). 
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The  modeled  tilt  response,  B(k),  matches  the  data  fairly  well  at  k  <  0.2  cpm 
(Figure  E3);  the  coherence  is  the  same  as  for  Hr{k).  Weakening  the  tail  force  term, 
Qt i  to  0.1  m-1  s_1  significantly  reduces  |£{&)|,  demonstrating  the  dominance  of 
forcing  on  the  tilt  response  (there  is  no  such  sensitivity  to  any  other  parameter).  For 
the  free  response,  approximated  using  Qn  =  Qt  =  0.01  m-1  s-1,  |B(fc)|  retains  the 
shape  of  the  weakened  Qj-  =  0.1  m_1  s_1  version,  but  is  at  a  fifth  its  level.  Past 
0.2-0.3  cpm,  the  model  is  inadequate  in  reproducing  the  observed  tilt  response. 

The  motion  and  orientation  of  the  vehicle  account  for  the  difference  between 
the  ACM  measurements,  ?Vm,  and  the  ocean  velocity,  uwat.  Our  model  infers  one 
portion — the  combination  of  tilt  position  and  oscillation — from  the  accelerometer 
measurements  and  manufactures  the  other  portion — the  point-mass  motion — based 
on  the  modeled  horizontal  dynamics.  Their  relative  contributions  are  compared  in 


Figure  E3  Comparison  of  model  tilt  response  (heavy  line)  with  that  inferred  from 
PATCH  EX  accelerometer  measurements  (dash-dot  line);  agreement  is  good  out  to 
0.2  cpm  and  fair  out  to  0.3  cpm.  Light  line  is  B{k)  for  weakened  tail  forcing,  com¬ 
puted  with  Qj  =  0.1  in-1  s-1  and  with  remaining  parameters  at  their  optimal  values. 
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Figure  E2  using  the  transfer  function  (Cr  —  Cm)jy\.  from  umsp  to  —BjWlliSV0+Lvd0ldt. 
Point-mass  motion  dominates  tilt  effects  for  k  <  0.1  cpin.  Between  0.16  and  0.33  cpm, 
they  are  closely  matched  both  in  amplitude  and  phase  (within  ±30°).  Tilt  oscillation 
rises  at  0.4  cpm  and  becomes  the  dominant  mode  by  0.6  cpm.  These  observations 
have  implications  on  the  sensitivity  of  the  results  (cWttl)  to  the  values  of  Bj  and  Lv\ 
there  is  little  effect  from  changes  in  Bj  since  they  are  felt  mainly  at  k  <  0.15  cpm; 
effects  from  Lv  are  fairly  even  at  all  k  >  0.1  cpm,  so  the  shape  of  the  resulting  ocean 
shear  spectrum  doesn’t  change  significantly  from  minor  alterations  of  this  parameter. 

4.  Dominance  of  Tilt  in  Accelerometer  Data 

To  check  our  assumption  that  we  cau  reduce  (10)  to  (13)  by  ignoring  all  but 
the  tilt  term  in  the  accelerometer  measurements,  we  form  the  ratio  l)(k)/B(k)  (Fig 
ure  E4).  For  the  data,  our  use  of  (13)  implies  a  constant  amplitude  of  0.8  m  s-2  and 
a  phase  of  zero.  The  deviation  of  the  modeled  response  at  k  <  0.1  cpm  is  due  to 
the  lateral  acceleration,  dvUUiV/dt.  This  is  of  no  concern  because  the  ultimate  results, 


Figure  E4.  Transfer  function  from  vehicle  tilt  to  measured  acceleration,  i.e.,  D(k)/ B{k) 
Heavy  line  is  model  response  for  optimal  large-tail  parameters.  Light  line  corresponds  to 
the  simplified  6  —  ameas/y  used  for  tin*  data  and  has  a  constant  amplitude  of  9.X  m  s“2 
and  zero  phase.  Bashed  line  shows  response  with  tail  forcing  weakened  by  setting 
Q,  -  0.1  nr1  s-‘  . 
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oceanic  velocities  or  shear  spectra,  are  insensitive  to  tilt  effects  at  large  scales.  The 
fluctuations  past  G.(i  cpm  are  artifacts  of  the  modeled  point-mass  response  (for 
which  is  unreliable  and  dominated  by  tilt  oscillation  in  this  band.  The  downward 
trend  at  k  >  0.1  cpm  comes  from  the  centripetal  term,  Lu  (PO/dt1,  which  implies  that 
we  may  he  slightly  underestimating  0  at  small  scales.  Without  the  strong  tail  force, 
the  relation  between  tilt  and  measured  acceleration  becomes  more  complicated,  as 
illustrated  by  the  dashed  line  computed  with  Q- j-  =  0.1  m-1  s-1. 


The  magnetometer  detects  vehicle  tilt  but,  unlike  the  accelerometers,  is  insen¬ 
sitive  to  the  dynamics  of  vehicle  motion.  Tilt  angitis  computed  with  (Al)  from  the 
magnetometer  data  are  consistent  with  those  computed  with  (13)  from  the  accelerom¬ 
eter  data.  They  are  almost  perfectly  coherent  and  in  phase  up  to  0.5  cpm,  where  the 
magnetometer  data  begin  to  get  noisy  (Figure  105).  This  close  agreement,  coupled 
with  the  results  of  the  model  analysis,  supports  our  use  of  (13)  for  determining  till 
angle's. 
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Figure  E5.  Comparison  of  northward  till  computed  from  accelerometer  data  with 
that  computed  from  magnetometer  data  for  five  PATCIllOX  drops. 


E8  TR  9414 


UNIVERSITY  OF  WASHINGTON  •  APPLIED  PHYSICS  LABORATORY. 


5.  Summary 

Expressions  for  model  transfer  functions  are  useful  for  understanding  the  behav¬ 
ior  of  the  existing  instrument  or  for  planning  modifications  or  new  construction.  We 
were  able  to  find  parameters  to  model  MSP’s  tilt  response.  Although  the  model  re¬ 
produces  the  observed  tilt  only  at  larger  scales,  it  is  useful  for  judging  the  sensitivity 
of  the  response  at  other  scales  as  well.  We  fine  tuned  Bj  and  Lv  with  the  transfer 
function  //,„  from  accelerometer  to  ACM  measurements.  We  showed  that  tilt  effects 
are  important  in  interpreting  ACM  data  only  at  scales  less  than  10  m  and  that  the 
final  results  are  weakly  sensitive  to  minor  variations  in  Bj  and  Lv.  We  demonstrated 
the  success  of  computing  tilt  directly  from  accelerometer  data,  especially  for  those 
scales  where  the  ACM  measurements  are  affected.  This  simplified  approach  is  pos¬ 
sible  for  MSP  largely  because  free  oscillation  is  damped  or  dominated  by  the  forced 
tilt  response. 
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